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This work deals with the evolution of a FORTRAN 
simulation written by Naval Research Laboratory which is 
used to evaluate the effectiveness of deploying ship- 
launched infrared decoys to counter the anti-ship infrared 
seeking missile threat. Although the model (referred to as 
the Stochastic Infrared Engagement Model — SIREM) possesses 
extensive analytical capability and flexibility, refinements 
are desired to more accurately emulate atmospheric effects 
on the acquisition process. Methods are derived herein to 
calculate atmospheric tr ansmi ttance as a function of range 
using an accurate, LOWTRAN— based empirical formula. 

Basic seeker d i scr i mi nat i on techniques are addressed 
which may be incorporated into SIREM or other simulations 
for future mi ssi 1 e— ver sus-decoy evaluations; and some ideas 
are presented which may prove valuable in decoy enhancement 
to subjugate the d i scr i mi nat i ng seeker. 
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I. INTRODUCTION 



With the advent of the missile era has come a revolution 
in naval warfare methods- The problem of contending with 
"smart" weapons has given a whole new meaning to both 
offense and defense alike- These missiles may employ a 
variety of acquisition hardware. Some will be passive, some 
active, and some both. 

Particular problems arise when missiles use passive 
acquisition methods- Passive sensors radiate nothing and 
are therefore invisible to conventional electronic support 
measures (ESM) . Furthermore, their small size allows them 
to reflect only meager amounts of radar energy, which 
further reduces reaction and coun termeasur e employment time 
of the ship. The countermeasures may be in the form of 
onboard or off board weapons systems or deception devices. 

The cost versus effectiveness of onboard close-in weapon 
systems has been a driving factor in the Navy's desire to 
analyze the employment of offboard expendable seduction and 
distraction devices as a relatively inexpensive alternative. 

Here we will restrict our studies to the passive 
infrared seeking missile problem, although concurrent 
analysis is being conducted to contend with the active RF 
missile seekers. 
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To reduce radically the cost of development of off board ♦ 
device character i st i cs and deployment tactics, a computer 
simulation has been constructed to assist in the analysis 
process- This modelling program is appropr i atel y named the 
"Stochastic Infrared Engagement Model - SIREM". It 
originated at Naval Research Laboratory in Washington, D.C., 
and was transported to the Naval Postgraduate School in 
Monterey, CA to assist in the analysis of the use of 
infrared decoys as a counter to the passive seeker problem. 

It is designed to imitate the typical cruise-missile versus 
ship engagement scenario and to establish measures of 
effectiveness (MOEs) for both 1) decoy char acter i st i cs , and 
2) deployment tactics. Therein lies the basis for 
continuing verification and upgrades to the model so that it 
remains capable of accurate evaluation of the hardware at 
hand. 

Most real— world systems are very complex. In order to 
model any system accurately, all factors which affect it 
must be represented in as accurate a fashion as possible. 
Weaknesses in repr esentat i on should be identified and 
documented so that overconfidences in the output do not 
occur. Once documented, the weak areas of the model may be 
studied and reworked to increase the accuracy of the model. 

This paper deals with the evolution of SIREM. By 
simulating each component of the engagement scenario in a 
realistic manner, the optimum decoy deployment schemes can 
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be obtained, combined with their RF counterparts, and 
algorithmically automated tor deployment via computer in 
order to reduce chances of operator control errors. It can 
easily be seen here that it is not -financially or 
logistically -feasible to evaluate all scenarios in -field 
tests using real ships and missiles, since the cost would be 
prohibitive. 

Extensive conversion and review o-f the SIREM model was 
conducted to assist in its enhancement. Two specific areas 
o-f weakness were identi-fied which require some attention. 

The -first area o-f concern was a lack o-f atmospheric 
attenuation on the emitted IR signatures o-f the ship and/or 
decoys. The second regards the study of seeker head 
di scr i mi nat i on techniques and how to implement them based on 
data available. The emphasis in this paper revolves around 
-finding a simple technique -for calculating atmospheric 
transmittance throughout the path of an incoming infrared 
seeking missile which will accurately a-ffect its probability 
of acquisition in the simulation environment. Two methods 
currently used in obtaining tr ansmi ttance in a stochastic 
model such as SIREM involve either look-up tables which are 
either incomplete or inflexible, or i ncor por at i on of the 
LQWTRAN program as a subroutine. The second method is 
preferred since it is very flexible to observed atmospheric 
conditions; however, it requires an extension of 9000 lines 
of FORTRAN code along with the added input— output coding. 
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which serves only to calculate one value of atmospheric 
transmittance -for each incremental range throughout the 
missile flight path. It was felt by the author that this 
was excessively wasteful of valuable computer time and 
memory requirements considering the model's current size of 
approx i matel y 15,000 lines of code. 

The intent was to devise a simpler scheme for finding as 
accurately as possible what the atmospheric tr ansmi tt ance 
would be for a given set of atmospheric conditions that 
would be easily incorporated into the incremental flight 
path of a missile and more accurately affect the seeker head 
acqui si t ion probabi 1 i ty. 

To more fully understand the principle of atmospheric 
transmittance and its effects, one chapter is dedicated to 
explaining the basic principles of atmospheric 
transmittance. Once a good understanding of the principles 
of atmospheric tr ansmi t tance has been gained by the reader, 
an introduction to the LOWTRAN program is performed, 
presenting its capabilities for prediction of transmittance. 
This will be followed by a simplistic approach for modeling 
atmospheric tr ansmi tt ance with very short computer coding 
schemes. 

The closing chapters will be dedicated to more detailed 
explanation of the SIREM modelling program, followed by an 
introduction to possible di scr i mi nat i on techniques Since 
seeker d i scr 1 mi nat i on adds to the problems surrouding the 
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use of off board countermeasures, the subject of enhancing 
current decoy designs by taking advantage of certain natural 
phenomenon is addressed in the concluding Chapter- 

Future enhancements to SIREM to incorporate both 
transmi tt ance and di scr i mi nat i on techniques are the object 
of this work and concurrent work by Naval Research 
personnel- These enhancements are intended to more 
accurately reflect the acquisition probability of the 
seekers being modelled so that refinements to the design and 
deployment of offboard coun ter measures are possible. 
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I I - ATMOSPHERIC JEN VI RONM ENT 



This chapter is dedicated to reviewing the principles o-f 
atmospheric propagation o-f m-frared radiation, the causes oT 
its variability, and tools available to predict how much is 
transmitted at various ranges and atmospheric conditions. 

The atmosphere is comprised o-f numerous gasses — 
nitrogen, oxygen, water vapor, carbon dioxide, methane, 
carbon monoxide, and ozone to name the major components. 
Additionally there are numerous suspended particles 
distributed in various densities throughout the earth's 
atmosphere. The major constituents o-f the atmosphere are 
nitrogen, comprising 78"/., and oxygen, which occupies 207. o-f 
the total volume. Water vapor makes up 17., while carbon 
dioxide constitutes only 0.047. o-f the atmosphere. The 
transmission of infrared radiation through the atmosphere 
depends on the meteorol ogi cal conditions, and hence varies 
with weather conditions and altitude. Only an approximate 
result can be gained through the theoretical study o-f the 
principles o-f atmospheric attenuation since no definitive 
conclusions can be made. The base o-f knowledge regarding 
atmospheric attenuation and the associated -formulae comes 
principally -from empirical means or approximations. 
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The primary factor in absorption in the infrared 
spectrum is water vapor. The wide variations in 
concentration and distribution of water vapor in the 
atmosphere are the primary reasons for the variance in 
absorption levels from one meteorological profile to the 
next. Although carbon dioxide comprises only 4/100 percent 
volume of the atmosphere, it is second in importance as an 
attenuator of infrared radiation, as well as being more 
uniformly distributed than water vapor. 

There are three major factors which contribute to the 
attenuation of infrared radiation: Ex t i n e t i on , or the loss 

of energy resulting from the interaction of the beam with 
the various absorbing or scattering constituents contained 
in the atmosphere; Ref r act i on , or the bending of light rays 
due to the refractive index gradients of the atmosphere, and 
S c int i ll at i on , or the distortion of the optical beam due to 
small scale turbulence. Refraction will not be discussed 
heavily here since it will be apparent later that it may 
serve only to extend the optical horizon of infrared seeking 
missiles, but have little affect on the simulation discussed 
herein. Extinction represents the combined effects of the 
scattering of radiation into or out of the beam due to 
aerosols and air molecules present in the beam path, and 
absorption caused by the direct transfer of energy from 
radiation generated by the source to the vibrational or 
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rotational kinetic energy of the molecules and aerosols. 
This is illustrated in the following diagram. 
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Figure 2-1. Causes of extinction 

The attenuation or extinction of electromagnetic 
radiation is described by the Lambert-Beer-Bouguer Law 



I = I 0 



(Eg. 2.1) 



where I is the attenuated radiation, I Q is the source 
radiation, p is the attenuation coefficient, and R is the 
path length. The extinction is caused by absorption and 
scattering by molecules and aerosols. They contribute 
linearly to the the total extinction as follows: 



P — Pn a + Pns + P AA + Pa b (Eq. 2.2) 

where p n A = molecular absorption by water vapor , ozone, etc. 
Phb— molecular (Rayleigh) scattering 
P*a= aerosol absorption (dry particles neglected) 

Pa 8= aerosol (liie) scattering 

The importance of each of these effects depends on the 
wavelength of the electrooptic energy. The wavelength 
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dependencies exist because of the relative sizes of the 
responsible constituents (molecule or aerosol ) and the 
energy states of the molecules (molecular absorption bands). 
The importance of the principal component of attenuation for 
the primary wave bands is given here. 

TABLE 2.1. PRIMARY ATTENUATORS IN IR SPECTRUM 

Wavelength Region Attenuation Coefficients 

( In micrometers ) In order of i mportance 
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The aerosol relationship used to describe extinction is 
principally governed by the size (aerosol radius) of the 
particle and the wavelength of the incident energy as: 

1) Rayleigh Scattering (where r<<\ — molecule effects). 

2) Mi e Scattering (where r # X approx i matel y ) . 

3) Geometric Scattering (where r >>\ - large particles). 

An illustration of the scattering description is as shown 
here. 
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The attenuation scattering cross— sect 1 on tor one particle of . 
size r and wavelength \ is the integral of the scattering in 
all directions: 

S = f I (0) dfi (Eq. 2.3) 

J Sol 1 d Anglo 

Mie scattering is generally concentrated in the forward 
direction as follows. 



I (0) 




Note that a minimum occurs near the ±100° points and that a 
secondary peak occurs in the backward direction (back- 
scatter ) . 

Rayleigh scattering is symetrical as shown. 



180 ° 




$=G° 



Figure 2—4. Rayleigh Scattering 

The scattering efficiency Q of a particle is given by the 
cross-section S (above) divided by the area of the particle. 



Q = S 

nr 2 (Eq. 2.4) 

For Rayleigh scattering Q is on the order of X -4 , and based 
on the wavelength span present from the sun's radiation, is 
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the principal factor in the sky's blue appearance. 



Q is 



higher at the longer wavelengths. For geometric scattering 
Q # 2.0. Only particles with radii between about 0.1 and 
1.0 micrometers are effective in scatter^ rig light in the 
visible range. These particles are too small to rain out 
due to their low collection efficiency. This is why rain 
alone does not reduce haze and increase visibility. A 
charge in the vertical temperature profile (such as the 
breaking up of an inversion 1 ayer by a frontal passage) does 
increase horizontal visibility through upward diffusion of 
the aerosols. 

Because the scattering efficiency Q is a function of 
aerosol size and wavelength, scattering itself is a function 
of wavelength. 



Here (X is the Angstrom coefficient, and varies generally 
from 0.5 to 2.0. This value may be as high as 4.0 for 
Rayleigh scattering. The scattering coefficients for 
various aerosol types are shown below. 



8 £ X-* 



(Eg. 2.5) 
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Figure 2—5 Scattering Coefficient 
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With this diagram we see that the wavelength dependence of 
scattering explains why the sky is blue, since the shorter 
wavelengths are scattered so much more. Likewise heavy haz 
can cause the red sunset effect since the longer wavelength 
are not scattered out of the path. Note also how the cloud 
fog, and haze coefficients carry over into the infrared 
region. We can expect these to be major factors in 
tr ansmi tt ance reduction there. 

Scattering is also indirectly a function of relative 
humidity (RH) because as RH increases, the water molecules 
collect on a dry particle until enough of them are on the 
particle to cause it to go into a solution and form a 
droplet. This is called deliquescence. Since a solution 
droplet can grow in size much faster than a particle, the 
scattering increases rapidly with RH due to an aerosol 
rapidly changing its size distribution through growth of 
droplets. Typical scattering coefficient variations are 
shown below. These are presented as a ratio of the 
scattering coefficient for aerosol s to the scattering 
coefficients for aerosols at 207. relative humidity as a 
function of relative humidity. 






(RH-20X) 




Figure 2—6 Scattering Coefficient Variation 



Models of the aerosol distribution in the marine 



environment can be very complex because of the amount and 
size of the particles present in the marine environment. 
These are functions of particle generation at the sea 
surface and the vertical transport of the particles once 
they are suspended in the atmosphere. 

Wind speed dependence on aerosol distribution in the 
maritime model is similar to that of continental models 
until the wind speed nears 14 knots. At 14 knots the 
aerosol concentration has increased by a factor of about 1.5 
due to aerosol generation by the surface whitecaps which 
form at this wind speed. As the wind speed increases, the 
concentrat i on may rise to a value of 10 to 100 times the 
calm atmospheric value near the surface, but decrease in 
concentrat i on with altitude up to appr ox i mat el y 5 
ki 1 ometer s. 

It can now be seen that the atmospheric conditions have 
a dramatic effect upon the extinction and absorption at 
various wavelengths. In this paper the interest lies 
primarily in the infrared region since this is where our 
model (SIREM) is designed to simulate target acquisitions. 
Hence it is good practice to become familiar with those 
factors which affect the systems being modeled. 

I ncorpor at i on of "real-world" occurrences is a vital part of 
systems modeling. Those models which use random (non— 
deterministic) methods to emulate a patterned parameter can 
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only be termed less than valid- They would serve useful for 
nothing more than gaming exercises since the answers could 
not be relied upon. 

Since the SIREM model is designed to obtain measures of 
effectiveness of IR decoy deployment versus an IR seeking 
missile, those measures may only be obtained through 
realistic i ncorporat i on of the factors that affect the 
engagement. The primary intent of the following chapters is 
to define the scenario, present the principles of detection, 
and show how the LOWTRAN program can be used to develop a 
very compact and accurate method of i ncor porat i ng 
atmospheric absorption into the acquisition process- 
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III. ENGAGEMENT 5CENAR I O 



As with any detection method, whether RF , visible, or 
infrared, it is assumed that the target is either emitting 
or reflecting energy at some power density in the frequency 
band of the detector as in the scenario shown below. 




From extended ranges this emission or reflection can be 
assumed to appear as appoint source, however at closer 
ranges the physical size of the target must be expanded into 
a distributed source emitter, especially considering the 
typical optical magni f i cat i on performed by the seeker. 

Since this document deals specifically with infrared 
emission and tr ansmi ssi on , the case of target reflection 
will henceforth be disregarded. 

The unit of measure given to the power emitted from the 
target is the radiant intensity (J), commonly stated in 
watts per steradian (solid angle). The target, assumed to 
be a point source, emits a source level radiation in 
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specific spectral bands. The radiant flux density of any 
wavelength across the continuum has a natural reduction in 
value based on the inverse square law with range from the 
ship (target) to the missile (receiving platform). 

The detector in the seeker head is composed of a 
material which electrically (via photon energy) sensitive to 
incoming radiation across some defined frequency bandwidth 
inherent to the composition. This bandwidth is preferrably 
fairly narrow. If it is allowed to become excessively wide, 
the si gnal — to— noi se ratio is reduced an unacceptable level 
of false alarm detections as discussed by Hudson CRef. 2:ch. 
12D. For this reason most detector systems have external 
spectral filters to narrow the bandwidth and eliminate these 
mi sdetect i ons. 

Based on the composition of the detector material each 
detector will have an optimum electrically active wavelength 
to which it is sensitive. As the incoming energy shifts in 
frequency to either side of the optimum, the electrical 
activity in the detector is reduced, creating some 
distribution which reflects the electrical activity, or 
molecular transition density as a function of wavelength. 
This relationship should be known in order to reflect any 
spikes in the sensitivity curve where certain energy bands 
may be particularly absorptive of incident energy. The 
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profile shown in Figure 3-2 is one which could be used in 
the typical simulation environment based on the above 
pr i nci p 1 es. 

Detector 
Radiation 
Detect i vi ty 
(Rel at i ve) 

2.0 3.0 4.0 X.. 

Figure 3-2. Detect i vi ty over Bandwidth 

The actual electrical bandwidth of a detector has three 
primary dependencies: first, the material from which it is 

constructed; second, the temperature at which it is 
operated; third, the physical area of the detector surface. 
The detectivity (D~) is the convenient measure categorizing 
a detector's operation within these parameters. It 
normalizes the "sensitivity" of various detectors to an 
electrical bandwidth of 1 Hz , and an area of 1 square 
centimeter, and is a convenient means for comparing the 
equivalent electrical response between detectors in circuits 
of different bandwidth. Shown in Figure 3—3 is a typical 
Lead Sulfide detector's D~ char acter l st i cs taken from Hudson 
[Ref. 2:pp. 3651. 
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Lead Sulfide (PbS) 



Type: Photoconduclive. 

I lenient Horn 0 0| x 0 01 to 25 x 25 mm. 
sizes Rectanguljr elements Ironi 0 0 1 
available: mm vside f oi an«iys. minimum 
element spacing is 0.0 1 mm. 
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Figure 3-3. Detectivity versus Wavelength 



The value of D~ is dependent upon the temperature, 
wavelength and bias current applied. The low temperatures 
are required, as well as bias current, to maximize the 
system's si gnal — to— noi se ratio. 

The object of our simulation is to characterize the 
search and acquisition process as closely to actual as 
possible. To do this we must model the detection process 
(seeker), atmospheric effects, flight character i st i cs of the 
missile, and target (ship or decoy) signatures as closely as 
possible. Taking all factors into account, the maximum 
range of the detector must be calculated. 

The spectral irradiance from a target (assumed distant 
in the field of view (FOV) is 

H v = Ta( \) in Watts/cm 2 (Eq. 3.1) 

R 2 
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where T A (X) is the tr ansmi ttance of the path, and the 
subject of -following chapters. R is the distance from the 
seeker to the target. 

The spectral radiant power into the detector is 

P x = HxAoTo(X) in Watts (Eq. 3.2) 

where Ao is the area of the entrance aperature of the 
optics, and T 0 <X) is the tr ansmi ttance of the whole sensor 
system. T P (X> is the product of the transfer functions for 
the detector and all external filter functions such as 
protective windows, lenses, reticles, and maskings. 

The signal voltage from the detector is 

V 8 = Pxft<X> (Eq. 3.3) 

where &(X) is the spectral responsivity of the detector. 

Thus far the equations represent an i nf i n i tesi mal 
spectral interval centered about some wavelength X- Since 
there is a finite bandwidth to an actual system, this should 
be taken into account as shown in Hudson CRef. 2:ch. 133. 

Remember, however, that we are assuming "worst case" in the 
simulation process. So, by selecting the optimum detectivity 
wavelength of the particular detector, we have effectively 
maximized the range of detection for our seeker. This 
assumption may only be made if the seeker system being 
modelled has a relatively narrow spectral response due to 



24 



external -filtering for a better signal-to-noise ratio. This, 
is virtually always the case from design inception since the 



false alarm problem is critical here. 

Further development of the equations above is performed 
by Hudson and leads to the Optical Trade-Off Equation: 



R = CJT a D 1 ' 2 |TD 0 (NA) 



1 1 D 0 (NA) T 0 1 1/2 C D~ I 1/2 [ 1 l l/2 

[ 2 I [ (WAf ) 1/2 <V 8 /V N ) I 



\ 1 / 



/ \ 



J (Eq. 3.4) 



1 — Target irradiance and atmospheric tr ansmi ttance 

2 — Optics parameters 

3 — Detector character i st i cs 

4 - System parameters and signal processing 



where R 
D 0 
NA 

Af 

v 8 /v N 



detection range. 

diameter of entrance of aperture optics, 
numerical aperture CRef. 2:ch. 51. 

i nstantaneous field of view of sensor, 
frequency bandwidth in Hz. 
signal to noise ratio LRef. 2:pp. 419D. 



By breaking the equation into sections as shown above we 
see in the first term that reduction of target signature 
(ship) and enhancement of the decoy signature inband would 
shift the acquisition probability more toward the decoy. In 
the engagement process nothing can be done about the 
atmospheric transmittance, although it does have dramatic 
effects on the acquisition range. Knowing the parameters of 
the second, third and fourth term facilitates modelling 
various missile/seeker systems. 
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At ranges beyond the sensitivity of the seeker head, 
i.e. where incident energy is insufficient, detection 
probability would be zero- In addition to the inverse 
square law power loss, the atmospheric losses previously 
discussed must be taken into account in order to accurately 
determine detection probability- Since the detectivity of 
the seeker is independent of the atmospheric effects on the 
emitted signature, the resultant probability of detection 
would be the product of probability of detection assuming no 
atmospheric effects, and the atmospheric tr ansmi ttance which 
varies from 0 to 1. 

Few models if any actually emulate the signal 
integrating principle upon which most detectors depend for 
threshold detection. The scanning technique, whether serial 
or parallel, performs the temporal integral of the incoming 
energy across the "window function" which defines the 
acceptance bandwidth of the seeker. The seeker systems 
probability of detection curve-ver sus-wavel ength is 
proportional to the detector's detect i vi ty-ver sus-wavel ength 
curve. 

The present version of the SIREM model emulates a 
scanner type seeker with limited field of view typical of 
the ASMs previously discussed. This model is presently 
being upgraded to account for atmospheric losses over range. 
It will do so by assuming the "worst case" posture. By 
this, it is meant that at the optimum detector wavelength. 
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the model presently calculates probability o-f detection -for 
incremental missile movements- This value is calculated as 
a "worst case" value tor the ship based on the detectivity 

0- f the seeker and the maneuvering capability o-f the missile. 
By multiplying the probability o-f detection (assuming no 
losses) at each range increment by the correspondi ng 
transmittance value at that range, the target acquisition 
probability is more accurately represented- This is the 
case whether the acquired target is a decoy or the ship. 

As an example assume that the transmittance has been 
calculated -for a given atmospheric pro-file -for the 
engagement scenario- The transmittance is then known at 
each increment o-f range. As the missile approaches its 
target the atmospheric transmittance increases approximately 
exponentially- I-f the probability o-f detection is 
calculated -for the seeker system in each incremental 
movement per -formed by the simulation, the engagement 
scenario probability o-f detection -for the missile would be 
the product o-f transmittance and original seeker acquisition 
probability generated by its own separate technique- Hence 

1- f the original P* were 0-8, and the transmi ttance were 0-6, 
the total system probability o-f detection would then be 
0 . 48 . 
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IV. 



ATMOSPHERIC MODELLING 



The previous chapters presented the principles of 
propogation of infrared radiation throughout the atmosphere 
and the factors which affect the acquisition process. This 
chapter will deal with an empirical approach to modelling 
the atmospheric tr ansmi tt ance using very accurate single 
equations once the conditions of visibility, range and wind 
speed are given. 

Atmospheric tr ansmi ttance losses for infrared radiation, 
not unlike radio frequency transmission losses, exhibit the 
character i st i c 1/ (Range) 2 loss factor as shown below. 




An actual curve varies widely in its rate of transmission 
loss with range depending heavily upon the particular 
infrared wavelength in question, the amount and size of 
particulate matter in the air, the altitude(s), and weather 
conditions as discussed previously. The LOWTRAN program 
CRef. 13 was developed to a large extent from empirical 
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observations and is used to predict atmospheric 
transmittance -for any given set o-f conditions. The desired 
goal in the continuing development of the SIREM model was to 
find a simple technique for i ncorporat i ng the calculation of 
atmospheric tr ansmi ttance into the incremental range 
movements of a simulated incoming missile and its associated 
infrared seeker head in order to more accurately portray the 
probability of acquisition of the seeker. 

A. LOWTRAN 

The LOWTRAN code calculates both atmospheric 
tr ansmi ttance and radiance throughout the infrared spectrum 
(0.25 to 28.5 mi crometers) . It uses a single parameter band 
model for molecular absorption, and includes the effects of 
continuum absorption, molecular scattering, and aerosol 
absorption. Refraction and earth curvature are included in 
the calculation for slant paths. 

Five seasonal models and the 1962 U.S. Standard 
Atmosphere are provided as inputs to the LOWTRAN program 
using typical altitude, pressure, temperature, water vapor 
density and ozone density profiles for each. The five other 
atmospheric profiles consist of (1) Tropical (15° N) , (2) 

Midlatitude Summer (45° N-July), (3) Midlatitude Winter (45° 
N-January) , (4) Subarctic Summer (60° N-July), and (5) 

Subarctic Winter (60° N-January) . With these profiles the 
user may select the general atmospheric conditions 
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applicable to the particular simulation. Additionally the 
user may input radiosonde data if desired to make 
calculations for a specific profile. Particular information 
regarding mixing ratios and assumptions made by the LOWTRAN 
program can be obtained from the LOWTRAN Manual CRef. ID and 
the associated texts listed therein. 

The LOWTRAN code includes three boundary 1 ayer aerosol 
types: Rural, Urban , and Maritime. Due to the nature of the 

SIREM/ASM model, only the latter was used since the 
composition and distribution of aerosols of oceanic origin 
are peculiar to the IR/ASM problem. These aerosols are 
largely sea-salt particles in the lower boundary layer 
caused by evaporation of sea-spray droplets. Together with 
a background aerosol of more or less pronounced continental 
char acter i st i cs they form a fairly uniform maritime aerosol 
which is repr esent at i ve of the lower 2 to 3 kilometers of 
the atmosphere over the oceans. 

These aerosols should not be confused with the heavy 
direct sea— spray aerosols found in the lower 10 to 20 meters 
above the ocean surface which is strongly dependent on 
recent wind velocities and sea surface turbulence. The 
extinction and absorption coefficients are calculated as a 
function of the relative humidity and based heavily upon the 
refractive index for the given altitude and weather 
condi t i ons. 
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The ability to calculate slant-path transmittances is 
incorporated utilizing a data base of seasonal and latitude 
dependencies of vertical aerosol distributions- This allows 
transmittance calculations of "hi gh-f 1 yer n missile attack 
profiles or aircraft infrared missile engagements by 
changing the input parameters to the LOWTRAN program to 
match those of the simulation requi rements - Iterating this 
calculation throughout the slant path range values yields 
the data points necessary to model the reduction in 
tr ansmi ttance over the entire engagement— range scenario with 
a single equation. The horizontal path used by the SIREM 
model for simulating cruise missiles is assumed to be a 
constant pressure path where earth curvature and refraction 
effects are negligible, or taken into account through the 
use of the 4/3 earth radius model. An alternative flat 
earth may be assumed for short range situations. 

The LOWTRAN code calculates transmi t tance as a function 
of "absorber density" for the path, the pressure, and 
temperature, for the particular wavelength band chosen. It 
utilizes both empirical laboratory data and available 
molecular line constant data in performing its cal cul at i ons. 
The output for a typical LOWTRAN tr ansmi ttance run is shown 
in Figure 4—2. The input to the program for spectral range 
is in the form of wavenumbers (inverse centimeters). In 
this example the range is 900 to 1145 cm’ 1 at increments of 5 
cm’ 1 required by the LOWTRAN program. The lesser number of 
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INTEGRATED ABSORPTION FROM 900 fO 11*.S Cn- 1 * 1 1 3 . 6 0 A VE R A GE I RA h SM T T A nC E * .S363 



wavenumbers is always the -first input into the LOWTRAN 
program, while the higher wavenumber is the second input. 
This corresponds to the hi gh-wavel ength to low-wavelength 
output in column 2 of Figure 4-2. The conversion to 
wavenumbers may be made using the following formulas. 



X = 10~ 2 / ( Wavenumbers) or 
Wavenumbers = 10“ 2 /X 



(Eg. 4.1) 



Here the wavenumbers are in inverse centimeter values as 
used by the LOWTRAN program. Tr ansmi ttance values for each 
of the contributing absorbers (H 2 0, C0 2 , Ozone, N 2 , Molecular 
Scattering, Aerosol, and Nitric Acid) are then calculated 
for each wavenumber (wavelength in column 2) in row format. 
The total tr ansmi t tance for each wavenumber is then 
calculated as the product of the tr ansmi ttances of each of 
the absorbing elements and is listed in column 3. The 
average tr ansmi ttance listed at the bottom of Figure 4-2 is 
simply the average of the tr ansmi ttances in column 3, and 
reflects the "integrated" overall tr ansmi ttance across, in 
this case, the 8.7 to 11.1 micrometer bands for the chosen 
input parameters (model, visibility, range, etc.). 

Since the bandwidth of most IR detectors is limited, 
either by its own material limitations or external filters, 
the input wavenumber range will seldom be as wide as the 
previous example. The 3.5 micrometer band was selected for 
the development of this paper and consists of wavenumber 
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inputs of 2850 (lower) to 2860 (upper) cm" 1 , and correlates 
to 3.496-3.509 micrometers in 5 cm" 1 intervals. 

In addition to the output shown in Figure 4—2 the 
program generates a large quantity of "header" information 
for each run, which reprints the inputs to the program. The 
unfamiliar user should obtain hardcopy of early runs to see 
what the output consists of. Once comfortable with the 
location of various values of the printout, use of the video 
screen is more appropriate and less wasteful of paper. The 
objective of each run is the "average transmittance" for the 
given input parameters, which will later be used in the 
curve fitting process. 

One can easily see that in order to curve fit the 
tr ansmi ttance as a function of range with each other 
parameter held constant, enough runs through the LOWTRAN 
program must be conducted at "selected ranges" to obtain a 
smooth curve fit. This means that for every model, range, 
wind condition, etc. the LOWTRAN program must be run in 
order to obtain the appropriate transmittance value. This 
would be a monumental task for the computer, and without 
stripping the PRINT statements from the LOWTRAN code, could 
generate a massive output file. The objective here is to 
show that a minimum number of selected ranges (7 or 8) can 
be used for this process to decrease computer time spent. 

The task becomes much easier if a particular set of 
conditions is in mind, i.e. a Midlatitude Summer model with 
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average wind conditions o-f 5 m/s and visibility o-f 30 
kilometers tor example- This circumstance could be modelled 
by running a "stacked input" -facility o-f the LOWTRAN program 
containing the ranges desired- A test run can be performed 
to -find the approximate maximum range where the 
transmi ttance value -falls below some operator desired 
threshold; e.g- T =0 - 0025 . 

Curve -fitting atmospheric transmittance as a -function of 
range requires that all other parameters (i.e. wind, 
visibility, etc-) remain constant- Additionally, enough 
points on the range scale must be taken to provide the 
necessary accuracy across the span of range of interest. 

The missile engagement scenario requires that transmittance 
be a factor in the acquisition process from horizon to 
target , thus the span in range should be from appr ox i matel y 
20 kilometers to zero- After extensive empirical analysis 
the values in Table 4-1 were selected to meet the above 
criteria. 

Up to this point it was assumed that all parameters 
(including wind and visibility) were constant- The desired 
goal is to analyze the individual effects of wind and 
visibility on transmi ttance and develop an empirical 
equation which generates a good approx i mat i on to the 
transmittance value for any wind and visibility condition 
over the scenario range, similar to that shown in Figure 4— 
1- The transmittance takes on a negative exponential 
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characteri stic whose decay constant varies as a function of 
wind and visibility. The idea is to empirically establish 
the decay constant as a function of both wind and 
visibility, the two principal parameters which must remain 
variable to the model user to real i st i cal 1 y model any 
engagement scenario. Intuitively (and verified by LOWTRAN) 
atmospheric tr ansmi t t ance decays much more rapidly if either 
visibility decreases or wind increases. To establish how 
the decay constant changes with wind and/or visibility, 
enough values for both wind and visibility must be analyzed 
to curve fit the decay constant dependence to either. The 
values of wind and visibility in Table 4.1 were empirically 
selected and adequately cover the conditions for most 
scenarios. Then, for each combination of wind and 
visibility, transmi ttance values over the range span are 
obtained using LOWTRAN. This builds a 3— d i mensi onal table 
which is tha basis for curve fitting and establishing the 
decay constant as a function of three variables; l-e. range, 
visibility, and wind. 

TABLE 4.1. RECOMMENDED CURVE FIT DATA POINTS 
Range in Km Wind m / s Vi si bi 1 i ty l n Km 



0.2 


0 . 1 


40 


1.0 


5.0 


30 


2.0 


10.0 


20 


4.0 


15.0 


15 


8.0 


20.0 


10 


12. 0 






18.0 






21.0 
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The method used -for data reduction is an extremely simple 
technique implemented -for mi crocomputer spreadsheets. The 
method employs a scheme -for performing Weighted Least 
Squares curve -fit techniques to the data, and while 
calculating the coefficients o-f the curve -fit, the goodness 
o-f -fit o-f the curve to the data may be viewed graphically by 
the modeller -for approval or correction as necessary. Some 
examples o-f the process will be given later. 

The LOWTRAN program serves as an excellent tool -for 
calculating transmi ttance over a wide variety o-f parameters, 
however, the program consists o-f 9000 lines o-f FORTRAN code. 
Most simulations are already quite large, and the addition 
of this code size as a subroutine is undesi reabl e , 
considering that each increment in range of the missile 
flight path must vector to LOWTRAN as a subroutine in order 
to obtain an updated single transmi ttance value. 

The alternative method presented here for obtaining the 
tr ansmi ttance value desired is by utilizing LOWTRAN external 
to the simulation and curve fitting the output data points 
as a function of range to a single equation, given the 
visibility and wind conditions of the particular engagement. 
This method requires extensive use of the LOWTRAN program to 
obtain the data points necessary to perform the curve fit. 
Once completed though, the two coefficients of the 
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approximated curve are the only data necessary to write an 
equation of the form: 

T = 0 Exp(-tfR) (Eq. 4.2) 

The coefficients of this equation ( tf and «) are all that is 
necessary to allow calculation of transmittance as the 
missile moves incrementally through the stochastic 
simulation. When an increment inrange occurs, a vector to 
the transmittance calculating subroutine (or function) with 
this equation updates the transmittance affect on the 
probability of acquisition- As previously discussed, the 
decay constant («) is actually a function of wind and 
visibility. Later in this chapter accomodations will be 
made for this. 

The SIREM model specifically calculates through 
simulation, the effectiveness of deployment of Torch/PIP 
Torch type infrared decoys against infrared seeking anti- 
ship missiles (ASMs) . It is assumed from this point on that 
the incoming missiles utilize a "low altitude" attack 
profile typical of the cruise missile variety which prevents 
early detection by the target ship. Assume for the time 
being that it flies at an altitude of approximately 60 feet 
(20 meters from the LOWTRAN example modelled herein). This 
would constitute a horizontal transmittance path when using 
the LOWTRAN program to calculate transmittance throughout 
the mi ssi 1 e-to-target flight path. Some other assumptions 
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utilized in the empirical derivations in this paper which 
were used as inputs to the LOWTRAN program are; 

a) The Marine Atmospheric Boundary Layer model was the 
only one utilized. 

b) A land haze -Factor o-f 3 was used (de-Fault). 

c) The wavelength modelled here was 3.5 micrometers. 

d) Each o-F the 6 atmospheric pro-files was calculated. 

e) De-Fault parameters inherent to LOWTRAN were used, 
varying only the visibility conditions, range, and 
wind conditions 

-f ) There was no prec i p i tat i on . 

Each of these parameters may be changed if desired to 
meet an alternative modelling requirement. For example, if 
it is known that a "new" IR missile threat has evolved which 
utilizes a "high— flyer" attack profile and a 3.2 micrometer 
optimized detector, the modeller could then choose the 
slant— path option in the LOWTRAN program selected at the 
necessary wavelength. Knowing only the flight profile, 
detector f r equency ( i es) and atmospheric conditions would 
then allow full use of the empirical tr ansmi ttance modelling 
technique presented here without insertion of 9000 lines of 
LOWTRAN code. 

Another advantage to this approach other than its 
simplicity and accuracy is that, if the simulation is 
written in a language other than FORTRAN, this method can 
still be easily adapted due to the short code length 
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required using the standard exponential function. The 
equational form is easily translated as will be seen. 

Let's analyze an example scenario. Suppose we are faced 
with our shi p-versus-mi ssi 1 e engagement described 
previously, and wish to calculate atmospheric transmi tt ance 
at various ranges in order to approximate its affect on the 
acquisition process. First we must assume some 
meteor ol og i cal conditions appropriate to our example for 
inputs to the LQWTRAN program. We will select here the 
par ameter s: 

1) Midlatitude Summer profile. 

2) Visibility is 30 kilometers. 

3) A maritime environment. 

4) A land haze factor of 3 (may vary from 1 to 10, with 
10 being bad haze as around Los Angeles). 

5) Wind is virtually nill. (0.1 m/s) and has been so 
for the past 24 hours. 

6) The range from the seeker to target is 4 kilometers. 

7) Assume a 3. 5 micrometer band (2850—2860 cm' 1 ) 

8) The path is horizontal at 20 meters altitude. 

We may now insert out data into the L0WTRAN program in the 
format listed in the L0WTRAN Manual CRef. 11 . 

In doing so we would find that at our 4 kilometer range, 
the tr ansmi ttance for this wavelength would be 0.6516 (See 
Appendix A: Table A. 2. a). The value may differ only 

slightly between computers. This means that the received 
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energy density at a range of 4 kilometers is only 657. as 
high as that for a perfectly transmitting medium (vacuum). 
We can now duplicate our efforts for the 4 kilometer range, 
but now with a wind velocity of 10 meters per second, 
keeping visibility the same. Doing this we obtain a 
tr ansmi ttance value of 0.4537. In the diagram below we see 
that increasing the wind, or decreasing visibility, has a 
dramatic effect on the decay constant associated with the 
tr ansmi ttance. 



Atmospher i c 
Tr ansmi ttance 



Fi gure 




Range (Km) 



4-2. Wind Effects on Transmittance. 



Note the approximately exponential fall off with range. 
This is the curve we wish to approximate as: 



T(R,V,W) = tf-EXP(0C (V) -f (W) -R-6) (Eq. 4.3) 

from the data points obtained from the LOWTRAN program. 
Here T, R, V, and W represent transmittance, range, 
visibility, and wind respectively. There is correlation 
between wind and visibility, and the effects are accounted 
for in the factor 6. 
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In order to model any particular atmospheric profiles, 
the data points for the 3-dimensional table may be obtained 
by adapting the following overlay algorithm as the main 
program, using LOWTRAN as a subroutine. 



FOR Atmospher i c_Prof i 1 e - 1 to 6 
FOR Wind_Speed = 0 to 20 Step 5 
Wind = 0.01 + FLOAT ( Wi nd__Speed ) 
FOR 1 = 1 to Nr__of _Range_Val ues 
READ Range^Value 
OPEN Input_F i 1 e 
WRITE Input_File 
CALL LOWTRAN 
WRITE Output 
NEXT I 

NEXT Visibility 
NEXT Wind_Speed 
NEXT Atmospher i c_Pr of i 1 e 



* See Note 1 

* See Note 2 



* See Note 3 



Note Is Prior to entering the looping structure the user 
must assign values to the other parameters used by the 
LOWTRAN program (e.g. wavenumbers, haze factor, etc.). 

Note 2: The actual wind speed may not be zero. If so 

LOWTRAN defaults to 4.1 m/s. Hence the following statement, 
Note 3: The input file must be formatted as specified in 

the LOWTRAN Manual CRef. 13. 
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Again it should be emphasized that the source code of 
the LOWTRAN program should be commented out to prevent 
printing of the massive administrative data. Only the 
"Average Transmittance" over the sufficiently narrow band is 
desired. The output may be directed to a printer or file. 
The latter method is preferred if input to the curve fitting 
process is desired in a direct fashion. In order to keep 
track of which data value is associated with which range, 
visibility, and wind value, it may be wise to write those 
associated values in the same output row. 

Continuing with our example we will assume for now that 
only one value for wind, and likewise visibility are 
calculated. The tr ansmi ttance would then become only a 
function of range since all other parameters remain 
constant. The desired goal is to select range values which 
best facilitate the curve-fitting process- Using the range 
values from Table 4.1 we can now perform an example of the 
curve fitting process, while presenting the spreadsheet 
technique used for doing Weighted Least Squares. 
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Examp le 



1 . 



Making the assumption that all parameters remain 
constant as before, with only range changing, a curve fit of 
the following LOWTRAN data is performed for the Midlatitude 
Summer profile with wind at 5 meters per second, and 
visibility of 30 kilometers. 

TABLE 4. 2a. LOWTRAN DATA FOR CURVE FITTING 



RANGE 


TRANSMITTANCE 


WEIGHTING FACTOR 


0.2 


0.9272 


8 


1 . 0 


0. 7824 


1 


2.0 


0. 6588 


2 


4.0 


0. 4828 


21 


8. 0 


O. 2723 


6 


12. 0 


0. 1588 


1 


18.0 


0.0727 


1 


21.0 


0. 0498 


2 



The Weighted Least Squares procedure is conducted in 
two parts- To effect the weighting, implementing multiple 
occurrence of the same data points is shown in table 4-3- 
Following in figure 4- 4a and 4- 4b is the actual Least 
Squares algorithm implemented in spreadsheet format. The 
output desired is the values for the A-Coef f i c i ent (O') and 
B-Coef f icient (a) for equation 4-2- In this example the 
first coefficient 0=0-875824, and the second, <X=-0. 140899 - 
Using equation 4.2 yields the table which follows at the top 
of the following page- 
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TABLE 4.2b. DATA POINTS EXTRACTED FROM FITTED CURVE 



range 


TRANSMITTANCE 


0.2 


0.8515 


1.0 


0.7607 


2. 0 


0.6607 


4.0 


0.4985 


8.0 


0.2837 


12.0 


0. 1615 


18.0 


0 . 0693 


21.0 


0.0454 



Comparing table 4.2a with 4.2b (i.e. the LOWTRAN 

data versus the exponential Curve-Fitted data) we notice 
comparable transmi ttance values, especially at ranges 
greater than one kilometer. Due to the nature of the 
engagement scenario, we may assume transmi t tance values at 
ranges less than one kilometer or so to be equal to 1.0, 
since its effect on acquisition at these ranges is minimal. 
Above 1 kilometer the error is generally one to two percent 
at most, and hence, very useable for "probabi 1 i ty " 
modelling. Within the following tables the Least Squares 
process in implemented into spreadsheet format. Although in 
this case Lotus Development Corporation's 1—2—3 (TM) was 
used, the technique is similar if using commercially 
available packages. The data is arranged in row-versus- 
column format with the Least Squares algorithm implemented 
in Tables 4.4a and 4.4b. 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 



A B 



Model: Midi at Sum 

Visibility: 30 

Varied param: Visibility 



Range (km) Transmittance 



0.2 


0.9272 


0.2 


0.9272 


0.2 


0.9272 


0.2 


0.9272 


0.2 


0.9272 


0.2 


0.9272 


0.2 


0.9272 


0.2 


0.9272 


1 


0.7824 


2 


0.6588 


2 


0.6588 


4 


0.4828 


4 


0.4828 


4 


0.4828 


4 


0. 4828 


4 


0.4828 


4 


0.4828 


4 


0.4828 


4 


0.4828 


4 


0.4828 


4 


0. 4828 


4 


0.4828 


4 


0.4828 


4 


0.4828 


4 


O. 4828 


4 


0.4828 


4 


0.4828 


4 


0. 4828 


4 


O. 4828 


4 


O. 4828 


4 


0.4828 


4 


0.4828 


8 


0. 2723 


8 


0.2723 


8 


0.2723 


8 


0.2723 


8 


O. 2723 


8 


0. 2723 


12 


0. 1588 


18 


0.0727 


21 


0.0498 


21 


0.0498 



TABLE 4.3. 

TRANSMITTANCE AND RANGE 
DATA POINTS USED 



NOTE: 

Using the same data point 
multiple times at any point 
minimizes error about that 
point in the curve tit 
process. This emulates the 
weighting matrix commonly 
used in performing Wieghted 
Least Squares curve fitting. 
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TABLE 4.4a. CALCULATIONS 


USED IN LEAST 


SQUARED FITTING 




C 


: D 


: E 


s F 


1 


SUM of R, 


This section -finds the 


coefficients 


2 




A and B 


-for the equation: 


3 


210.6 


T = tf*Exp(tt*R) at 


Visibility =30 km 


5 


R, 2 


Ln (T.) 


CLn (T 4 ) D 2 


Ri*Ln (T» ) 


6 


0.04 


-0.0756 


0.0057 


-0.0151 


7 


0.04 


-0.0756 


0.0057 


-0.0151 


8 


0.04 


-0.0756 


0.0057 


-0.0151 


9 


0.04 


-0.0756 


0.0057 


-0.0151 


10 


0.04 


-0.0756 


0.0057 


-0.0151 


11 


0.04 


-0.0756 


0.0057 


-0.0151 


12 


0.04 


-0.0756 


0.0057 


-0.0151 


13 


0.04 


-0.0756 


0.0057 


-0.0151 


14 


1 


-0.2454 


0.0602 


-0.2454 


15 


4 


-0.4173 


0. 1742 


-0.8347 


16 


4 


-0.4173 


0. 1742 


-0.8347 


17 


16 


-0.7282 


0.5302 


-2.9126 


18 


16 


-0.7282 


0. 5302 


-2.9126 


19 


16 


-0. 7282 


0. 5302 


-2.9126 


20 


16 


-O. 7282 


0.5302 


-2.9126 


21 


16 


-0.7282 


0.5302 


-2.9126 


22 


16 


-0.7282 


0.5302 


-2.9126 


23 


16 


-0.7282 


0.5302 


-2.9126 


24 


16 


-0.7282 


0.5302 


-2.9126 


25 


16 


-0.7282 


0.5302 


-2.9126 


26 


16 


-0. 7282 


0. 5302 


-2.9126 


27 


16 


-0.7282 


0.5302 


-2.9126 


28 


16 


-0. 7282 


0. 5302 


-2.9126 


29 


16 


-0. 7282 


0.5302 


-2.9126 


30 


16 


-0. 7282 


0.5302 


-2.9126 


31 


16 


-0.7282 


0.5302 


-2.9126 


32 


16 


-0. 7282 


0.5302 


-2.9126 


33 


16 


-0. 7282 


0.5302 


-2.9126 


34 


16 


-0. 7282 


0.5302 


-2.9126 


35 


16 


-0. 7282 


0. 5302 


-2.9126 


36 


16 


-O. 7282 


0.5302 


-2.9126 


37 


16 


-0.7282 


0.5302 


-2.9126 


38 


64 


-1.3009 


1.6922 


-10. 4068 


39 


64 


-1.3009 


1 . 6922 


-10.4068 


40 


64 


-1.3009 


1.6922 


-10. 4068 


41 


64 


-1.3009 


1 . 6922 


-10.4068 


42 


64 


-1.3009 


1 . 6922 


-10.4068 


43 


64 


-1.3009 


1 . 6922 


-10.4068 


44 


144 


-1.8401 


3.3860 


-22.0813 


45 


324 


-2.6214 


6. 8718 


-47. 1855 


46 


441 


-2.9997 


8.9984 


-62.9945 


47 


441 


-2. 9997 


8.9984 


-62.9945 


48 


Sum (Ri ) 2 


Sum Ln ( T t ) 


Sum Squares 


Sum of Products 


49 


2079. 3200 


-35. 2421 


49-9966 


-320. 8972 
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TABLE 4 -4b. SPREADSHEET LEAST SQUARES CALCULATIONS 



G : H 



1 : Scratchpad 

2: -0.1409 

3: -0.1326 

4: 20.3153 

5: 20.4251 

6 : 

7: A— Coefficient 0.875824 
8: B-Coef f i c i ent -0.140899 
9: R-Squared 0.9946 



: Formula in cell at left: 

= Note: Here n=42 

> (42*F49— D49*C3> / ( 42*C49-C3^2 ) 

> (D49— G2*C3> /42 

> G2* (F49-C3*D49/42) 

> E49— D49*~2/42 

> Exp ( G3) 

> G2 

> G4/G5 



This satisfies the Least Squares process. Cells G2 through 
G5 represent respectively: 

Gl) CnE (R t 1 n (Ti ) ) - El n ( T* ) ER* D / CnER t 2 - (ER t ) 2 D 

G2) C El n (T t ) - ( Gl ) ( ER* ) 2 1 /n 



(Eq. 4. 4-4. 7) 

G3) G2 C ERi 1 n ( T t ) - ER t El n(T t )/n] 

G4 ) El n ( Ti ) 2 - ( El n (Ti ) ) 2 /n 

Cell C3 is the sum of the range values of column "A". 

C3) ER t (Eq. 4.8) 

The R— Squared value represents the "Goodness of Fit" 
of the function to the original data when using the 
approx i mat i ng equation. In this case the fit is a very one, 
but in cases of either extremely good or poor visibility 
(i.e. >40 or <5 kilometers) errors between L0WTRAN 

predictions and empirical results grow, and the fitted curve 
will have regions where error may be on the order of five to 
six percent. For modelling in the probabi 1 i st i c sense this 
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is still adequate, since the region of error is transient 
and errors of this magnitude have only small effects on the 
acquisition probability. Figure 4-3 is a plot of the above 
example using both LOWTRAN— obtained data points, and the 
correspondi ng exponential curve fit data points obtained by 
using the spreadsheet technique presented. We may now 
extend our process to calculate the atmospheric 
transmittance for any of the 6 LOWTRAN profiles for any wind 
or visibility conditions. 



Figure 4 — 3 



Comparison — LOWTRAN vs Exponential Fit 
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a. Extending the Curve Fit 

The 3-dimensional table (A. 2. a) of transmittance 
values has been obtained for, in this case, the Midlatitude 
Summer Profile- Thus far we have completed the curve- 
fitting process for one particular wind and visibility 
condition (i.e. visibility of 30 Km, wind speed 5 m/s as 
calculated in Table 4.4). The A and B coefficients (0 and 
a), which were derived as a result of Least —Squares fitting 
an exponential function, were used in equation 4.2 to 
emulate extinction over range for this wi nd /vi si bi 1 i ty 
combination to obtain Figure 4—3. There the accuracy of the 
exponential fit versus the original LOWTRAN data is visually 
apparent. Tables A.l.a through A.6.b in Appendix A list the 
values of tr ansmi ttance obtained from the LOWTRAN program 
and from the curve fitting process for each of the six 
atmospheric profiles available- The curve-fit coefficients 
for the other five atmospheric profiles are contained in 
Tables A. 7 and A. 8 in Appendix A. 

We now wish to extend this procedure to be able 
to emulate extinction over range for any wind/visibility 
combination. To do so we must first build a table of the A 
and B coefficients and examine how their values change from 
one wi nd/ vi si b l 1 i ty combination to another. Using the 
combinations from Table 4.1 the coefficients were obtained 
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for the entire atmospheric profile of our example. These 
are contained in Table 4.5. 

By examining Table 4.5 we see that the A 
coefficients hardly change throughout entire range of wind 
and visibility values. This coefficient may be assigned a 
constant value for the remainder of the extinction curve-fit 
process. This is true when modelling each of the 
atmospheric profiles individually; however, the A— 
coefficient does change from one profile to another. This 
can be seen in Table A. 7 in Appendix A. 



TABLE 4.5. EXPONENTIAL CURVE FIT COEFFICIENTS 
FOR "MIDLATITUDE SUMMER" PROFILE 



A— Coef f i c i ent Val ue 



Wi nd 


Vi s=40 


Vi s=30 


Vi s=20 


Vi s= 15 


Vi s=10 


0 


0.875788 


0. 875808 


0. 875754 


O. 875834 


0.875812 


5 


0.875747 


0.875824 


0.875699 


0.875477 


0.877011 


10 


0. 875820 


0. 875716 


O- 875650 


0 . 874656 


0.874052 


15 


0-375753 


0.875685 


0. 875366 


0. 87601 1 


0.872420 


20 


0.875777 


0 . 875885 


0.875886 


O. 877431 


0.859462 



B-Coef f i ci ent Val ue 



Wi nd 


Vi s=40 


Vi s=30 


Vi s=20 


Vi 3=15 


Vi s=10 


0 


-0. 064275 


-0. 065945 


- 0 . 0693 1 3 


-0.07271 1 


-0. 079494 


5 


-0. 118579 


-0. 140899 


-0. 185840 


-O. 231112 


-0.322254 


10 


-0. 129858 


— O. 156395 


-0.209959 


-0. 263693 


-0.371488 


1 5 


-0. 132919 


-0. 160663 


-0.216479 


-0. 273163 


-0.384626 


20 


-0. 134394 


-0. 162714 


-0.219772 


-0.277832 


-0.386820 



The B— coef f i c i ents do vary throughout Table 4.5. 
Through further examination, the change in its value is non- 
linear with wind and visibility, and as such, simple linear 
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i nterpol at i on of this table is not accurate. Nor is this 
technique desirable because it requires many more B— 
coefficient data points than those contained in Table 4.5 to 
be stored in memory to be even moderately accurate. The 
process of extending equation 4.2 to the more general form 
of equation 4.3 becomes a matter of finding out how the B— 
coefficients in Table 4.5 change as a function of wind and 
visibility (i.e. find the exponential decay constant of 
extinction in terms of its dependence on wind and 
visibility), and fit functions which will accomodate the 2— 
dimensional decay constant variation. 

Thus the curve fit process is continued, only 
now it must be performed on the B-coeff icients which vary 
with both wind and visibility. This will leave 'Range 4 as 
the only variable to input to equation 4.3 for calculating 
transmittance. From table 4.5 it can be seen that the B— 
coefficient increases negatively (decreases) as either 1) 
wind increases, or 2) visibility decreases, and is therefore 
a function of wind and visibility. Equation 4.3 stipulates 
that the decay constant for extinction over range consists 
of three parts. These are: 



a) 


« (V) 


— Due 


to 


changes in visibility. 


b) 


f <W) 


— Due 


to 


changes in wind, and 


c) 


<5 


— Any 


correlation between the first two 
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By also observing that the B-coef f icient is 
least negative where wind is minimum and visibility is 
maximum, we can now use the ratios of any other B— 
coefficients relative to the first and find their 
mul tiplicati ve increase over the original (baseline - 
maximum visibility, minimum wind) B-coef f i ci ent due to wind 
and visibility- This done by dividing each value in the B— 
coefficient table by the 'baseline' value in the upper left 
corner of the table as is done for Table 4-6- Table 4-6 
depicts how much greater (negatively) the decay constant is 
for any particular wind and visibility combination than it 
would be for its best transmittance (i-e. wind zero, 
visibility maximum)- The non-linearity of the B-coef f i ci ent 
with wind or visibility is much more evident in Table 4.6. 

We shall see in a few moments why this table is necessary. 

TABLE 4-6- RATIO OF B(Vis,Wind) TO B ( Vi s„ AX , Wi nd MIN ) 



Wind 


Vi 5 
40 


Vis 

30 


Vis 

20 


Vi s 
15 


Vis 

10 


0 


1.0 


1 . 02598 


1 . 07838 


1. 13125 


1 . 23678 


5 


1 . 84487 


2. 19213 


2.89133 


3 . 59567 


5.01367 


10 


2.02035 


2.43322 


3.26657 


4. 10257 


5.77966 


15 


2. 06797 


2. 49962 


3.36801 


4. 24991 


5.98407 


20 


2.09092 


2.53153 


3.41924 


4.32255 


6.01820 



The following steps are the procedures used to 
obtain the three functions of the exponential decay constant 
of equation 4-3- Individually these steps produce the 
functions a (V) , f (W) , and 6 which will allow the input of 
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any wind or visibility value -from 0-20 m/s or 40—5 Km 
respectively, and generate the approximated B-coe-ff icient of 
extinction for that condition. Keep in mind that this 
example only covers the 'Midlatitude Summer' profile. The 
other profiles are approximated in the same fashion, with 
the results included in Table 4.7. 

(1) Step One. As visibility decreases the 
B-coef f i ci ent becomes more negative, causing transmittance 
to decrease at a faster rate (more extinction). Testing of 
various functions was conducted to generate the B- 
coefficient as a function of visibility. It was assumed to 
start with that wind is zero so that the curve fit process 
uses the B-coef f i ci ents of the first row on Table 4.5. 
Graphically the data points appear as in Figure 4—4. 



Figure 4 — 4 
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Care must be taken in selecting a curve-tit 



method to approximate any of the three -functions «(V), -f (W) , 
or 6 since the B-coef-f icient data points are negative. This 
eliminates exponential or logarithmic -function Least-Squares 
techniques since the natural logarithms o-f the data points 
are necessary, and are undefined. 

To approximate « (V) , a cubic function was 
selected of the form: 

Row Ctt (Vis) Iz A, + A 2 -V + A 3 - V 2 + A 4 -V 3 (Eq. 4.9). 

This eliminated the negative data point problem and produced 
highly accurate results for the few data points on hand. 

The values found for the coefficients of the cubic are 
contained in Table 4.7. These values were used to generate 
the fitted curve in Figure 4—4, while the original B— 
coefficients were plotted as the single point values. 
Remember: fitting of this function is performed only on the 

first row of B-coef f icients in Table 4.5. 

(2) Step Two is to establish a function which 
reflects the change in the decay constant with wind: i.e. 

f (W) . Graphically the data points plot as shown in Figure 
4-5. At first the function appears as an exponential decay 
which would be easily modelled. But, since the data points 
are negative this method is not easily performed. This is 
where the concept of Table 4.6 is useful. Given that the B— 
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Figure 4—5 




Figure 4 — 6 



f(Wlnd) Data Points — Translated 
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coefficient for any visibility condition (with wind zero) 
can be generated by completing step one, then f (W) becomes 
merely a multiplier to that value generated. 

To find out what to multiply (X (V) by to get 
the decay constant for wind conditions which are not zero, 
we fit a function to the first col umn of Table 4-6. Now, 
all the data points are positive, and any function is fair 
game. Graphically, the newly translated function appears as 
in Figure 4-6. After unsuccessful (inaccurate) attempts to 
fit a hyperbolic tangent (tanh) function to this curve, a 
much more accurate one was discovered: 

Column Cf (Wi nd ) 3 : 1 + ED 1 - ( 1 -Ex p <D 2 - Wi nd ) ) 3 (Eq. 4.10) 

(3) Step Three . Figure 4-6 shows a family of 
five curves which plot each of the five columns of Table 
4-6. The curve with the least maximum represents the first 
column in Table 4.6. The others in increasing fashion 
represent decreasing visibility conditions. The correlation 
factor 6 is necessary to increase the extinction when wind 
and visibility are not zero and maximum respectively. In 
this case 6 acts as a multiplier also, driving the decay 
constant more negative when its value becomes greater than 
one. Fitting a function to points vertically spaced on 
Figure 4—6 is the key to finding the correlation coefficient 
6. A function which works extremely well in this case is: 
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6 = C, - V c * 



(Eq. 4.11) 



The same effect may be obtained by fitting 



equation 4.11 to the data points of one of the lower rows of 
Table 4.6 if. they have been normalized a second time in the 
horizonta l direction <i.e. divide each coefficient in one 
lower row of Table 4.6 by the coefficient in the first 
column of that row). 



process for each atmospheric profile, leaving eight constant 
coefficients which, when incorporated into equation 4.3, 
yield an admirable approx i mat i on of transmittance values at 
range values between one and twenty kilometers throughout 
the entire atmospheric profile. This procedure was used to 
establish the constants for each of the six LOWTRAN 
profiles. These constants are included here for reference. 



This procedure completes the curve fitting 



TABLE 4.7 



CURVE FIT COEFFICIENTS FOR ALL SIX LOWTRAN PROFILES 



TROPICAL 

PROFILE 



MIDLATITUDE 

SUMMER 



MIDLATITUDE 

WINTER 



A, ! -0. 139154 
A 2 ! 6. 2SQ7E-03 

A 3 I -2.2354E-04 
A 4 ! 2- 609 IE— 06 

C, ! 13.74105 

C 2 1 -0.716806 

D, I 0.895 

D 2 ! -0.30 

0 i 0.880 



-0. 1249852 
6. 296102E— 03 
-2. 207060E— 04 
2. 619482E— 06 
15. 18516 
-0. 743405 
1 . 10 
-0.27 
0.875 



-0.08956 
6.31 731 7E— 03 
—2. 247950E— 04 
2. 627775E— 06 
29.876 
-0.923857 
2. 46 
-0.295 
0.9383 
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TABLE 4.7 (Continued) 



SUBARCTIC 

SUMMER 



SUBARCTIC 

WINTER 



STANDARD 

ATMOSPHERE 



A, ! -O. 109255 
A 2 ! 6.2545E-03 

A, ! -2.2228E-04 

A, ! 2. 5953E-06 

C, ! 19.1121 

C 2 ! -0.80492 

D, 1 . 440 
D 2 I — O . \jO 

o : o. 



-0.07970 
6. 4225E-03 
—2. 2857E— 04 
2. 6723E-06 
37. 730 
-O. 98561 
3.980 
-0.2825 
0.9636 



-O. 091037 
5. 6255E-03 
—2. 0021E-04 
2. 3407E-06 
24. 130 
-O. 86706 
1.812 
-0.34 
0.9208 



Note here that eight of the nine constants 
in each atmospheric profile eventually constitute «, the 'B' 
or extinction coefficient. The next step is to implement 
and explain the computer coding process. This procedure is 
contained in appendix B. We can now turn our attention to 
the SIREM model and examine the acquisition process. 
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V. 



THE STOCHASTIC INFRARED ENGAGEMENT MODEL (SI REM ) 



As discussed in chapter one, the most inexpensive means 
of developing equipment char acter i st i cs and tactics for 
their use is through simulation. The use of computers to 
analyze effectiveness of systems or subsystems is becoming 
the predominant development method for this reason. The 
SIREM model was designed to investigate the ef f ect i veness of 
deploying decoys which emit an infrared signature against 
infrared seeking anti— ship missiles (ASMs) . 

By examining through simulation such parameters as burn 
time, deployment altitude, deployment angle, decoy drift and 
descent char act er i st i cs if airborne, and distance from the 
ship, optimum design requirements and deployment tactics can 
be developed. As with any software model, to maintain its 
validity it must be periodically reviewed and upgraded to 
remove any synthetic qualities and replace them with 
accurate repr esent at i ons . 

There are ongoing efforts to incorporate atmospheric 
losses, di scr i mi nat i ve seeker capab i 1 i t 1 es , and last but not 
least the multi— ship mul t i — mi ssi 1 e analysis capability. 
Coincident with the work being done on the SIREM model, 
extensive investigation is being conducted regarding optimum 
deployment character 1 st 1 cs for chaff, repeater decoys, and a 
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variety of other off board countermeasures in order to 
integrate the use of offboard devices into a multi— band 
defensive system- To fully understand the SIREM model 
requires extensive programming experience and/or simulation 
expertise- Likewise an operational knowledge of available 
devices is advisable; however, a general description of the 
program is in order. 

A- THE MODEL 

The SIREM model calculates probabi 1 i t i es associated with 
distraction of an IR homing missile by a pattern of decoys 
launched by the attacked ship. Rather then the Monte Carlo 
method of simulation, a stochastic approach was utilized in 
order to significantly reduce computation time- The program 
uses realistic input parameters pertinent to the scenario 
under study, and evaluates statistical information about the 
distraction effectiveness which may be output in graphical 
or numeric form. 

The user may easily modify any of the input parameters 
under the following general headings: 

1- Ship characteristics 

2- Missile char acter i st i cs 

3- Decoy character i st i cs 

4- Scenario and run char acter i st i cs 
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Although the original version transported from Naval 
Research Laboratory was designed to operate under the NLX 
NAMELIST facility of the Prime (TM) computer there, 
significant changes were required to adapt the input/output 
to non— NAMELIST type since it was not available at Naval 
Postgraduate School; and to modify the graphics subroutines 
and hardware to operate on the Digital VAX 11/7S0 (TM) and 
the DI— 3000 (TM) graphics software available there. 

Output from the program is available in the form of 
graphs or printed output. A summary of the available graphs 
is listed here: 

1. Probability of distraction conditioned on initial 
mi ssi 1 e-to-shi p di stance. 

2. Probability of distraction conditioned on initial 

azimuthal approach angle. 

3. Pattern of ship lock— on range as a function of angle. 

4. Probability of distraction conditioned on initial 
azimuthal approach angle and a single initial missile 
to ship distance. 

5. Angular regions where probability of distraction 
exceeds a threshold. 

6. Locations where threshold crossings of plot 5 occur. 

7. Probability of distraction conditioned on initial 
mi ssi 1 e— to— shi p distance and sector of initial 
azimuthal approach angle. 

8. Ship orientation, decoy deployment, and wind 
di recti on . 

9. Effects of atmospheric refraction on IR at various 
incidence angles. 
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Being stochastic by design, certain parameters may be 
considered random variables tor initialization. The initial 
position and orientation of the missile relative to the ship 
is a good case. It is assumed -for startup that the missile 
is pointed radially inward toward the ship, and hence the 
effectiveness computations would yield a lower bound for 
effectiveness of decoy distraction. This is the basis for 
the 'worst case' situation described in earlier chapters. 

The initial range and azimuth approach angle of the 
missile are discrete random variables. The initial range 
may be changed to appropriate kilometer values; while the 
initial range may also be modified, it is designed to run 
multiple missiles at the ship every 360/N T heta degrees. 

The seeker char acter i st i cs and the random variables 
assigned to emulate it can change from missile to missile. 

In this sense, randomness must be assigned to accomodate 
scanning directions, boresight error, field of view, and 
probability of acquisition on multiple targets (i.e. ship 
plus any decoys). The model performs these adequately, 
using significant but valid assumptions regarding the 
process. 

The maximum initial range of the missile to the target 
is R 0 at the time of launch of the first decoy. This is 
considered time t=0. It then takes R 0 /V„ seconds for the 
missile to reach the ship at the origin, with V M being the 
missile velocity. The sample time Tj interval determines 
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the maximum number of time samples taken N T ®. Ti is set 
small enough to obtain an accurate representat i on of the 
scenario. Since the missile is always aimed at ship before 
acquisition, a lower bound on distraction results, and if no 
decoy is acquired then the probability of distraction is 
zero for that given initial range and azimuth. 

The choice of time interval T t is dependent on the 
missile velocity, target velocity, and the scanner sweep 
time T s . Since the area of the scanner field of view on the 
sea surface is large relative to the distance either the 
missile or ship travels during the sweep time, T 8 may equal 
T 4 . As different conditions arise a necessity for time 
scaling may result. Developments have begun to incorporate 
different scanner parameters. 

The projection of field of view is incorporated to allow 
both curved earth and flat earth calculation. The ship and 
decoys are treated as point targets for infrared radiator 
purposes. The "hot spot" of the deployed decoy is placed at 
a known distance above the sea surface. Calculations are 
made to determine if the hot spot is within the field of 
view of the scanner. The next stage of development requires 
incorporating radiance and tr ansmi ttance calculation, and 
detector sensitivity parameters to determine if there is 
sufficient radiance from the ship or decoy to be detected by 
the seeker. Much of this development was covered in earlier 
chapters. 
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The model can also calculate any shadowing of decoys by 
the ship. This means that if the aspect of the ship comes 
between the missile and decoy, effectiveness of that decoy 
will inherently drop to zero. This can be done using either 
the flat or 4/3 earth model. In addition to shadowing by 
the ship, if excessive range exists, shadowing by the earth 
may occur. 

Once acquisition of a decoy occurs, the missile flies 
according to the aerodynamic parameters of the control 
surfaces and systems for whatever missile is being modelled. 
If a decoy is not in the field of view and the ship is 
acquired, the missile is assumed to hit the ship. If locked 
onto a decoy, the missile flies at the decoy until it passes 
it. At that point break-lock occurs. This condition may 
also occur if the control system forces the missile to 
maneuver so that the decoy is no longer in the field of 
view. And lastly, if the decoy burns out break— lock occurs. 

When break— lock occurs the scanner must re— enter the 
acquisition mode. Since the missile usually makes drastic 
maneuvers near break— lock, the re-acqui si t i on process is 
somewhat random, but if the ship is in the field of view it 
may be assumed to hit the ship. Missile flight continues at 
break-lock until all targets are passed. 

The flexibility of the program is excellent. Virtually 
any parameter is changeable either at input or at source 
code level. As many as five decoys may be deployed by the 
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ship (more if certain array dimensions are changed) in any 
combination of deployment angles, times, altitudes, arid 
active lifetime. The ship may change size and speed. The 
wind is a factor on decoy drift and missile flight. And the 
missile may change speed, altitude of approach, and control 
system or seeker parameters. The variations are endless, 
and analysis of the many combinations pertinent to real- 
world scenarios may be undertaken to make maximum effective 
use of offboard countermeasures. Additionally, by further 
extending not only SIREM but other similar models as well, 
evaluation of descr i mi nat i on techniques and optimization of 
decoy signatures become possible avenues of further 
development. With this in mind we turn our attention to the 
acquisition process and observed phenomena. 
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VI. I R SEEKER SYSTEM AND COUNTERMEASURES 



Since the principal interest here lies in passive 
seekers, it is convenient to separate these into three broad 
categories: (1) Scanning systems such as search, track, and 

imagers, (2) Staring systems which use non— mechani cal means 
to obtain spatial scene viewing, and (3) Measurement systems 
such as radiometers, spectr ometer s and interferometers. O-f 
interest here are the -first two categories. Scanning 
detection schemes are most common since staring systems 
usually require sophisticated mul t i — d i mensi onal array 
detectors and massive numbers o-f amplifiers to support the 
i ndi vi dual el ements. 

A. SCANNING AND STARING SYSTEMS 

Scanning systems sample the radiant intensity 
distribution in their -field o-f view using various 
horizontal / vertical techniques. The output -from the 
detector in this case will be a linear analog o-f the 
distribution or a simple indication o-f the presence o-f a 
target and its location. The object plane is usually 
assumed at in-finity and the detector responds to radiant 
levels within the instantaneous -field o-f view o-f its 
element (s). This gradient in radiance is translated to 
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electrical signals which are then filtered, amplified, and 
referenced to a baseline (missile axis) for inputs to the 
various azimuth and elevation control systems. For imaging 
systems the signal is translated from the object to its 
correspond i ng image plane by painting out the radiance 
distribution using optical deflection onto a film, or beam 
deflection to a cathode ray tube. Using internal processing 
techniques in missile seekers, the translation is used as 
input to the missile control surface system to center the 
highest radiant intensity in the center of the total field 
of view. The use of either spiral or raster techniques to 
accomplish slow mechanical scanning is most common. In some 
cases the relatively high speed motion at altitude provides 
the vertical trace capability for the system, while a 
spinning mirror provides the horizontal trace. Regardless 
of the scanning technique employed, major concern has 
recently been given to the i ncorporat i on of d i scr i mi nat i on 
processing into newer seeker systems. If (or, more 
appropriately, when ) this occurs, the effectiveness of using 
offboard decoys could be significantly reduced. Therefore, 
an examination of possible d i scr i mi nat i on methods is in 
order . 

In dealing with imaging systems there are two basic 
approaches to describe the sensitivity of the detector. The 
first approach is based on the detector figure of merit (D~) 
described earlier. It regards the detector as a noise 
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source of fixed level which is determined from the 
measurements of fixed conditions similar to those 
encountered by the detector when installed in its normal 
operating configuration inside the seeker. 

The other approach which is suitable only for BLIP 
(Background limited mode) detectors regards a random arrival 
rate of photons from the background as the ultimate source 
of detector noise. In such a case the detector noise power 
density is described using the laws of photon rate and 
detector quantum efficiency. When the number of incoming 
(target) photons is small in comparison to the number of 
background photons the approaches are equivalent, since BLIP 
D~ is determined by background photon rate. If the number 
of target photons is larger than the average number of 
background photons though, the noise due to the random 
arrival time of target photons dominates and D* then is 
meaningless. This case is the unusual one. 

Image forming scanners are not usually designed to 
respond to an absolute radiance level, i.e. not for a dc 
response- They are designed to follow changes in radiance 
as the field of view is scanned. This includes gradual 
changes in the radiance level cor r espondi ng to low 
frequencies. DC restoration can be used to control the 
intensity levels at the output. Without this form of 
automatic gain control (AGO periods of signal blackout or 
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saturation would occur due to high 1/f noise and amplifier 
drift at low frequencies. 

The signal of interest from the scene is the variation 
in irradiance H caused by the variation in temperature or 
emissivity of objects in the scene and may be expressed as: 



AH 



where W 
T 

€x 

W x 




in Watt s (Eq. 
cm 2 



o o 

= i nstantaneous field of view in steradians 
= temperature in degrees Kelvin 
= emissivity (wavelength dependent) 

= spectral radiant emittance in W/cm 2 



6 . 1 ) 



If most noise in the system is caused by the detector, 
it may be expressed in terms of "noise equivalent flux 
density" (NEFD) . This is the minimum signal irradiance 
capable of producing a peak signal to rms noise ratio (S p /Nr) 
of one. 

The ultimate goal is to find the actual si gnal —to— noi se 
ratio (S/N) generated by the scanner. This is given as: 

S = Variation of irradiance in scene (AH) i n Wa tts (Eq. 6.2) 
N NEFD 

The noise equivalent power (NEP) in watts of the system 
can be obtained from the figure of merit of the detector 
(D~) and is expressed as: 
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in Watts 



(Eq. 6.3) 



NEP = VA rf - Af 

D* 

where A<d = area of the detector in cm 2 
Af = system noise bandwidth 

D~ = detector figure of merit in cm • Hz l/2 /watt 

The signal power on the detector is equal to the optical 
collector area times the signal flux density- Therefore the 
system noise equivalent flux density is: 

NEFD = 4 VAc -A-f in Watts/cm 2 (Eq. 6.4) 

*D„ 2 D~€„ 

where D 0 = diameter of the aperture of the optical system 
€ 0 = optical efficiency 

The difference in temperature and emissivities (AT and 
At) of objects in the scene contribute to the differential 
radiant emittance (AH), but due to the small contribution of 
A€ for non— BLIP detector operation, its term in equation 6.1 
is commonly neglected- By substituting for the first 

integral term the abbreviated value for radiant emittance 
becomes then: 

AH = M-ATW, in Watts/cm 2 (Eq- 6-5) 

* 

The i nstantaneous field of view W is defined by the 
detector area A«d and the focal length fl as: 

W = A** in steradians (Eq. 6-6) 

(f 1 ) 2 
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By substituting equation 6.6 into equation 6.4 the noise 



equivalent flux density is reestab 1 i shed as: 

NEFD = 4F V* * Af in Watts/cm 2 (Eq. 6.7) 

iD 0 D*t 0 

where F is the f -number of the optics. 

The si gnal -to-noi se ratio generated by the scanning 
system is then (from equation 6.2): 

S = AT D 0 D ~€ 0 V¥ (Eq. 6.8) 

N 4F VET 

Some careful considerations are in order here since the 
AH actually seen by the sensor is modified by the 
tr ansmi ttance of the optical components and the atmosphere 
as discussed earlier. If the bandwidth considered is 
sufficiently narrow, tr ansmi ttance losses may be compensated 
for by using the average tr ansmi tt ance values for each while 
still producing accurate results. The narrow bandwidth was 
the premise for the development of tr ansmi t t ance modelling 
in Chapter IV. 

Sensitivity of an image forming scanner is expressed by 
its noise equivalent temperature (NET). This parameter 
refers to the temperature difference AT within a scene 
element (with a given emissivity and average temperature of 
e.g. 300°K) which can produce a change in the electrical 
signal level equal to the rms system noise. The noise in a 
well designed system is mostly detector noise. By 
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i ncorporat i ng the atmospheric transmittance into equation 
6.8 and solving tor AT, the NET is obtained as: 

NET = 4F l /At~ in "K (Eq. 6.9) 

WD 0 D~€ 0 T ft W J | 

This an expression of the sensitivity of the detector system 
as a whole. It makes the assumption that the scene area 
over which the temperature difference AT occurs is much 
larger than the field of view. 

From the design aspect a major emphasis must be placed 
on the value of D~ for a given detector since its value may 
vary as much as four orders of magnitude among common 
detector materials CRef 2s Table 7. ID. In addition, the D~ 
of a particular detector will vary considerably from the 
published data owing to such factors as manufacturing 
process, electrical bias, cooling, loading, aging and 
aerodynamic heating of the optical window. 

The electrical bandwidth (Af) of the system is 
determined by the scan rate. It is proportional to the 
reciprocal of the dwell time of an image point on the 
scanning detector surface. The amplifiers connected to the 
detector can be optimized to obtain a maximum signal-to— 
noise ratio, especially if the char acter i st i cs of the signal 
and noise are known. In this fashion a number of 
target /ob ject discrimination techniques may be devised. 
Optimum filter design is extremely difficult to perform 
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using only a -few components, though a simple band— pass 
-filter with a flat response and a low and high -Frequency 
cutoff is often adequate. The very high signal to noise 
ratios realized by using matched filter (correl at i on > 
techniques for processing are virtually impossible due to 
the passive nature of the systems, but by tweaking the high 
and low frequency roll-off a much improved background 
control and target definition may be accomplished with only 
a slight reduction of the S/N ratio. 

If a total field of view of Q steradians is scanned at a 
constant rate without redundancy and with an i nst antaneous 
field of view W steradians, then the scan is completed in a 
frame time of t B - The dwell time t d of a point on the 
detector is then t 8 -W/ft. As the total field of view is 
increased (either in azimuth, elevation or both) linearity 
in scan rate is difficult to achieve CRef 4:pp 7351. The 
off-center loss in dwell time must be compensated for by a 
scan efficiency factor € a . As the scan efficiency 
decreases, an increased electrical bandwidth is required. 

In this case the optimum bandwidth required is 
approx i mat el y : 

Af R = 5fl€ 9 in Hz (Eq. 6.10) 

4*t a 

Actually this refers to the upper cutoff frequency since it 
is not practical to design for dc response. The lower 
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cutoff frequency can usually be designed to handle down to a. 
few Hertz which allows reproduction of at least most of the 
low frequency features in the scene- If a linear (vertical 
usually) array of detectors is used, with N elements in the 
array, then Af R can be reduced proportionately. 

Equivalent noise bandwidth Af will usually be larger 
than required bandwidth Af R because of the inherent 1/f 
noise, amplifier noise, and R-C filter roll-off- In terms 
of scan rate, the noise equivalent temperature of a scan 
system is: 



NET = 



2 

ts 



D 0 D~W* N 



2 5€ s 
€ 0 T* 



in °K (Eq- 6-11) 



By breaking equation 6-11 into the above groups the 
correspondi ng components may be individually analyzed- The 
first group represents system performance requirements. The 
second group represents mostly the physical components 
available at design inception, and the third group 
represents efficiency parameters of materials used and the 
atmosphere. 

The principles of detectors and line scanning systems 
are covered extensively in References C2 j and C33- The 
basics presented thus far will serve as a conceptual 
approach to modelling seeker systems, and will serve as a 
baseline for devising di scr i mi nat i on techniques in the 
modelling process- Based on discussion thus far no 
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committment has been made with regards to any particular 
system, but either conceived or existing seekers which use 
the line scanner techniques operate under the previously 
ex plained principl es. 

1 - Filtering -for Pi scr i mi nat i on 

There are three -filtering methods which may be used 
in seeker system designs. These are: 1) Spectral, 2) 

Spatial, and 3) Temporal. Each method has its physical 
limitations -for implementation in di scr i mi nat i on as will be 
discussed shortly. In contrast, the i ncorpor at i on o-f 
di scr i mi nat i on into IR seekers leads to the necessity to 
analyze present coun ter measur es -for possible enhancement to 
subjugate any o-f these discrimination techniques. 

a. Spectral Filtering 

Since background radiant intensity is a primary 
consi der at i on to seeker design, methods to reduce it are a 
necessity. Background radiant intensity typically has peaks 
(maxima) at approx i matel y 0.5 and 10 micrometers caused by 
solar reflection and thermal sel f — emi ssi on respectively. 
Around 3 microns a region of minimum background radiance 
exists. Coincidentally, between the two maxima, most hot 
targets produce their peak radiant intensity. The simplest 
method of detection is to examine typical target spectral 
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character i st i cs, select a detector most appropriate for 
operation at the target's peak spectral points, and then use 
filtering to prevent the background radiant energy from 
entering the system. If the target of interest exhibits 
significant wavelength-dependent emissions, one method of 
di scri mi nat i on could be the use of notch (2-color spectral 
bandpass) filters to select specific wavebands. If the 
radiant energy in the passbands produced by the decoy is 
then less than that produced by the target (W/sr) , the decoy 
may be rendered ineffective in the seduction role. Thus 
much of the design of a decoy revolves around matching its 
spectral radiance to that of the ship it is being designed 
to protect. 

Most IR decoys are omni -di r ect i onal by nature, 
and exhibit the same spectral signature in all directions 
(di sregardi ng wind or sea-surface effects if floating). 
Target signatures are aspect dependent (intensity) as well 
as spectrally dependent. As the aspect changes, so does the 
intensity of certain spectral components- The inherent 
variability of target signatures in this respect makes it 
exceedingly difficult to establish di scr imination techniques 
since look-up tables and matched filter methods cannot be 
used . 
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b- Temporal Filtering 



Spatial -Frequencies produce time varying signals 
with a -Frequency content equal to the spatial -Frequency 
times the angular scan rate. In processing then, low 
-Frequency rejection and high -Frequency peaking can be used 
to eliminate background clutter. Only i-F the target and 
background signals can be adequately described can a spatial 
filter maximize the ratio of target signal to background 
signal. If mul t i -el ement detectors are used, simpler filter 
techniques are possible. The instantaneous field would be 
made as small as practical, and the upper frequency cutoff 
of the electrical filter is selected to correspond to the 
dwell time of a point on the detector. The low frequency 
cutoff is usually set at 1/3 to 1/4 of the upper cutoff 
frequency to block the background radiance without allowing 
saturation at the other end of the bandpass. This filtering 
technique is most often used to decrease the probability of 
false alarms through scan— to— scan comparisons and has little 
merit in discerning targets from decoys. 

c. Spatial Filtering 

A more promising approach to d i scr 1 mi nat i on by 
the seeker may be to process spatial character i st i cs of the 
scene. When the range from the seeker to the target is 
large, the objects in the scene will appear as point 
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sources- As this range decreases, however , large objects 
-form a distribution within the scene (both horizontal and 
vertical) while small objects retain their point source 
character i sties- By keeping the instantaneous -field o-f view 
small, integral (in spatial sense) processing could be used 
to discern large distributed objects from small ones- It 
then becomes the task of the decoy designer to not only make 
the decoy more intense spectrally than the ship, but to make 
the signature more widely distributed to simulate large 
targets. If the decoy is suspended in the air, for example 
by parachute (and somewhat less so if the decoy is 
waterborne) , the wind will tend to distribute the heat over 
a wider area; but variable wind conditions and short hang 
times of parachute-deployed decoys make this method of 
deployment less attractive than water— borne placement- Some 
of the effects of wind for signature expansion are still 
present with water— borne decoys. In addition to the torch 
emission above the sea surface there exists a reflective 
component off the sea surface referred to as the "Halo 
effect" ERef. 41. This effect was originally noticed by 
personnel from Naval Research Laboratory taking infrared 
measurements of the USS Ticonderoga's signature in both the 
3—5 and 8—14 micron bands. The halo effect may readily be 
seen in Figure 6-1. Upon investigation it was established 
that the sea foam caused by a white— capped sea surface has 
very high reflectance at low grazing angles in both infrared 
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Figure 6-1. 



Haloing caused by sea-surface 
reflection around USS Ticonderoga 
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bands. 



The white-caps may be caused by wind, ship movement , • 



and other surface turbulence. When viewed -from certain 
angles as in Figure 6-2, energy -from the source which was 
not originally emitted in the direction of the seeker 
(observer) is reflected from the sea surface toward the 
seeker, adding to the original signature radiant density. 

From the decoy design aspect this effect may be 
capitalized upon by i ncorporat i ng a means of disturbing 
(foaming) the sea surface over a wide area around a water- 
borne decoy beneath its burning source. By floating the 
decoy at the proper height above the foamed sea surface a 
very large reflective signature may be achieved in addition 
to the intense point source of the torch itself. 

Foaming may be accomplished in many ways, but 
primary interest must be in distributing over as large an 
area as possible. Since it is difficult to foam the sea 
surface over a wide area using a simple squirrel-cage fan or 
compressed air (too many cylinders) another prospect must be 
taken which considers weight and space restrictions of the 
decoy. Foaming of sea water is purely transient and very 
short in duration, so the above methods would probably be 
discarded upon testing. The alternate means suggested here 
is to use a small pneumatically driven suction pump with a 
mixing chamber to mix the incoming sea water with, for 
example, light— water or aqueous film-forming foam (AFFF) . 

The mixture may then be sprayed from the decoy by fan or 
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nozzles over a wide area of the sea surface to create a 
highly reflective 1 ayer (using sea surface turbulence for 
wide angular reflection) and corr espondi ngl y large 
signature- The components in this case could be very 
inexpensive- Additionally, this technique would work in all 
infrared bands and, because of the difficulty in devising 
passive-only di scr i mi nati on methods, may be the optimum 
offboard countermeasure to IR seekers- The methods 
discussed here have not been attempted, but these and many 
other approaches will need to be addressed in order to meet 
the IR threat with inexpensive offboard countermeasures. 
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TABLE A.l.a 



TRANSMITTANCE IN 3.5 UM BAND - TROPICAL PROFILE 
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Range 
in Kb 


*** 

40 


LOWTRAN DATA 

Visibility in 
30 20 


*** 

Ki lonet 
15 


ers 

10 


5 


0.20 


93. 17% 


93.14% 


93.07% 


93.01% 


92.88% 


92.50% 


1.00 


82. 05* 


81.91% 


81.63% 


81.35% 


80.80% 


79.10% 


2.00 


73.19% 


72.95% 


72.46% 


71.97% 


70.99% 


68.16% 


4.00 


60.33% 


59.92% 


59. 11% 


58.30% 


56.73% 


52.28% 


8.00 


43.64% 


43.06% 


41.90% 


40.77% 


38.59% 


32.77% 


12.00 


32.70% 


32.09% 


30.81% 


29.58% 


27.27% 


21.40% 


18.00 


22.12% 


21.46% 
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21.00 
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8.80% 
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7.28% 


8.00 


28.30% 


23.63% 


16.48% 


11.46% 


5.55% 


0.64% 


12.00 


17. 10% 


13.05% 


7.60% 


4.40% 


1.50% 


0 . 06% 


18.00 


8.30% 


5.56% 


2.47% 


1.09% 


0.21% 


0.03% 


21.00 


5.90% 


3.70% 


1.44% 


0.55% 


0.08% 


.00% 


0.20 


91.76% 


91.47% 


90.49% 


89.52% 


87.61% 


82.14% 


1.00 


76.84% 


74.83% 


70.92% 


67. 19% 


60.31% 


43.70% 


2.00 


64.21% 


60.88% 


54.70% 


49.10% 


39.56% 


20.78% 


4.00 


46.41% 


41.73% 


33.68% 


27.13% 


17.61% 


48.60% 


8.00 


25.83% 


20.88% 


13.60% 


8. 83% 


3.72% 


0.28% 


12.00 


14.90% 


10.84% 


5.70% 


2.98% 


0.82% 


0.02% 


18.00 


6. 80% 


4.21% 


1.61% 


0.61% 


0. 09% 


0.01% 


21.00 


4.67% 


2.67% 


0.87% 


0.28% 


0.03% 


. 00% 


0.20 


91.90% 


91.39% 


90.37% 


89.36% 


87.37% 


81.69% 


1.00 


76.61% 


74.51% 


70.46% 


66.59% 


59.49% 


42.51% 


2.00 


63.81% 


60.37% 


53.98% 


48.22% 


38.49% 


19.65% 


4.00 


45.85% 


41.03% 


32.80% 


26. 17% 


16.68% 


4.35% 


8.00 


25.21% 


20.19% 


12.91% 


8.22% 


3.33% 


0.23% 


12.00 


14.38% 


10.30% 


5.27% 


2.68% 


0.69% 


0 . 00% 


18. 00 


6.43% 


3.90% 


1.43% 


0.52% 


0.07% 


0.00% 


21.00 


4.38% 


2.45% 


0 . 76% 


0.23% 


0.02% 


0.00% 


0.20 


91.87% 


91.35% 


90.32% 


89.28% 


87.26% 


81.47% 


1.00 


76.50% 


74.36% 


70.24% 


66.31% 


59.11% 


41.96% 


2.00 


63.62% 


60. 12% 


53.65% 


47.82% 


38.00% 


19.15% 


4.00 


45.59% 


40.71% 


32.40% 


25.74% 


16.25% 


4. 13% 


8.00 


24.92% 


19.87% 


12.59% 


7.95% 


3. 17% 


0. 20% 


12.00 


14.08% 


10.03% 


5.08% 


2.55% 


0.64% 


0.01% 


18.00 


6.27% 


3.77% 


1.35% 


0.48% 


0 . 06% 


0.01% 


21.00 


4.23% 


2.33% 


0.71% 


0.21% 


0.02% 


.00% 
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TABLE A.l.b 



TRANSMITTANCE IN 3.5 UM BAND - TROPICAL PROFILE 

*** EMPIRICAL *** 



Range ! Visibility in Kilometers 



in Km 


40 


30 


20 


15 


10 


5 


Wind= 0 


0.2 


86.63% 


86.58% 


86.57% 


86.49% 


86.33% 


86. 


Wind= 0 


1.0 


81.37% 


81.13% 


81.08% 


80.72% 


79.94% 


78. 


Wind= 0 


2.0 


75.25% 


74.80% 


74.70% 


74.05% 


72.63% 


70. 


Wind= 0 


4.0 


64.34% 


63.59% 


63.41% 


62.31% 


59.94% 


56. 


Wind= 0 


8.0 


47.04% 


45.95% 


45.69% 


44.12% 


40.82% 


35. 


Wind= 0 


12.0 


34.40% 


33.20% 


32.93% 


31.24% 


27.80% 


22. 


Wind= 0 


18.0 


21.50% 


20.39% 


20.14% 


18.61% 


15.63% 


11. 


Wind= 0 


21.0 


17.00% 


15.98% 


15.75% 


14.36% 


11.72% 


8. 


Wind= 5 


0.2 


85.75% 


85.14% 


84.16% 


83.06% 


80.76% 


74. 


Wind= 5 


1.0 


77.30% 


74.60% 


70.42% 


65.94% 


57.28% 


38. 


Wind= 5 


2.0 


67.91% 


63.23% 


56.35% 


49.41% 


37.29% 


16. 


Wind= 5 


4.0 


52.40% 


45.44% 


36.08% 


27.74% 


15.80% 


3. 


Wind= 5 


8.0 


31.21% 


23.46% 


14.79% 


8.74% 


2.84% 


0. 


Wind= 5 


12.0 


18.58% 


12.11% 


6.07% 


2.76% 


0.51% 


0. 


Wind= 5 


18.0 


8.54% 


4.49% 


1.59% 


0.49% 


0.04% 


0. 


Wind= 5 


21.0 


5.79% 


2.74% 


0.82% 


0.21% 


0.01% 


0. 


Wind=10 


0.2 


85 . 55% 


84.88% 


83.82% 


82.63% 


80.13% 


73. 


Wind= 10 


1.0 


76.39% 


73.48% 


69.00% 


64.22% 


55.08% 


35. 


Wind=10 


2.0 


66.32% 


61.35% 


54.10% 


46.87% 


34.47% 


14. 


Wind= 10 


4.0 


49.98% 


42.77% 


33.26% 


24.96% 


13.50% 


2. 


Wind=10 


8.0 


28.38% 


20.79% 


12.57% 


7.08% 


2.07% 


0. 


Wind= 10 


12.0 


16.12% 


10.10% 


4.75% 


2.01% 


0.32% 


0. 


Wind= 10 


18.0 


6.90% 


3.42% 


1.10% 


0.30% 


0.02% 


0. 


Wind= 10 


21.0 


4.51% 


1.99% 


0.53% 


0. 12% 


0 . 00% 


0 


Wind= 15 


0.2 


85.50% 


84.82% 


83.74% 


82.53% 


79.99% 


73. 


Wind=15 


1.0 


76. 19% 


73.23% 


68.68% 


63.84% 


54.60% 


34. 


Wind=15 


2.0 


65.97% 


60.94% 


53.61% 


46.32% 


33.87% 


13. 


Wind= 15 


4.0 


49.45% 


42.19% 


32.66% 


24.38% 


13.04% 


2. 


Wind= 15 


8.0 


27.79% 


20.23% 


12. 12% 


6. 75% 


1.93% 


0. 


Wind= 15 


12.0 


15.61% 


9.70% 


4.50% 


1.87% 


0 . 29% 


0. 


Wind=15 


18.0 


6.58% 


3.22% 


1.02% 


0.27% 


0.02% 


0. 


Wind=15 


21.0 


4.27% 


1.86% 


0.48% 


0. 10% 


0.00% 


0. 


Wind=20 


0.2 


85.49% 


84.81% 


83.73% 


82.51% 


79.96% 


73. 


Wind=20 


1.0 


76. 15% 


73.17% 


68.61% 


63.76% 


54.49% 


34. 


Wind=20 


2.0 


65 . 89% 


60.84% 


53.50% 


46.19% 


33.74% 


13. 


Wind=20 


4.0 


49.33% 


42.07% 


32.52% 


24.25% 


12.94% 


2. 


Wind=20 


8.0 


27.66% 


20.11% 


12.02% 


6.68% 


1.90% 


0. 


Wind=20 


12.0 


15.50% 


9.61% 


4.44% 


1.84% 


0.28% 


0. 


Wind=20 


18.0 


6.51% 


3.18% 


1.00% 


0.27% 


0.02% 


0. 


Wind=20 


21.0 


4.22% 


1.83% 


0.47% 


0. 10% 


0.00% 


0. 



03% 

60% 

20 % 

01 % 

64% 

68 % 

52% 

21 % 

54% 

36% 

72% 

18 % 

11 % 

00 % 

00 % 

00 % 

41% 

56% 

37% 

35% 

06% 

00 % 

00 % 

00 % 

16% 

96% 

89% 

19% 

05% 

00 % 

00 % 

00 % 

11 % 

83% 

78% 

16% 

05% 

00 % 

00 % 

00 % 
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TABLB A. 2. a 



Wind 



Wind 



Wind 



Wind 



Wind 



TRANSMITTANCE IN 3.5 UM BAND - MIDLATITUDE SUMMER PROFILE 





*** 


LOWTRAN DATA 


*** 






Range . 




Visibility in 


Ki loieters 




in Kb 


40 


30 


20 


15 


10 


5 


0.20 


94.16* 


94.12* 


94.06* 


94.00* 


93.87* 


93.49* 


1.00 


84.47* 


84 . 33* 


84.05* 


83.76* 


83 . 20* 


81.53* 


2.00 


76.79* 


76.53* 


76.02* 


75.51* 


74.49* 


71.53* 


4.00 


65.59* 


65. 16* 


64.28* 


63.42* 


61.72* 


56.92* 


8.00 


50.26* 


49.59* 


48.27* 


46.98* 


44.50* 


37.84* 


12.00 


39.82* 


39.03* 


37.48* 


35.99* 


33.17* 


26.02* 


18.00 


28.86* 


28.01* 


26.36* 


24.79* 


21.94* 


15.24* 


21.00 


24.90* 


24. 04X 


22.40* 


20 . 86* 


18.09* 


11.82* 


0.20 


93.14* 


92.72* 


91.89* 


91.06* 


89.43* 


84.72* 


1.00 


80.01* 


78.24* 


74.80* 


71.48* 


65 . 29* 


49.83* 


2.00 


68.89* 


65 . 88* 


60.21* 


54.99* 


45.88* 


26.72* 


4.00 


52.78* 


48.28* 


40.33* 


33.64* 


23.41* 


7.94* 


8.00 


32.55* 


27.23* 


19.00* 


13.22* 


6.40* 


0.74* 


12.00 


20.75* 


15.88* 


9.26* 


5.37* 


1.81* 


0.07* 


18.00 


10.86* 


7.27* 


3.23* 


1.43* 


0.28* 


0.00* 


21.00 


7.96* 


4.98* 


1.94* 


0.75* 


0.11* 


0 . 00* 


0.20 


92 . 93* 


92.44* 


91.45* 


90.47* 


88.53* 


83.01* 


1.00 


79. 11* 


77.04* 


73.02* 


69. 18* 


62.10* 


45.00* 


2.00 


67.35* 


63.86* 


57.38* 


51.50* 


41.50* 


21.79* 


4.00 


50.46* 


45.37* 


36.62* 


29.50* 


19. 15* 


5.28* 


8.00 


29.74* 


24.05* 


15.66* 


10.16* 


4.29* 


0.33* 


12.00 


18. 13* 


13. 18* 


6.93* 


3.62* 


0.99* 


0.02* 


18.00 


8.87* 


5.49* 


2.09* 


0.79* 


0.11* 


0.00* 


21.00 


6.28* 


3.60* 


1.17* 


0.38* 


0.04* 


0.00* 


0.20 


92.87* 


92.36* 


91 . 33* 


90 . 30* 


88.29* 


82.55* 


1.00 


78.87* 


76.71* 


72.54* 


68.55* 


61.25* 


43.76* 


2.00 


66.94* 


63.32* 


56.62* 


50.58* 


40.37* 


20.60* 


4.00 


49.84* 


44.60* 


35.66* 


28.45* 


18. 13* 


4.72* 


8.00 


29.02* 


23.24* 


14.86* 


9.46* 


3.84* 


0.26* 


12.00 


17.47* 


12.52* 


6.40* 


3.25* 


0.84* 


0.01* 


18.00 


8.39* 


5.09* 


1.86* 


0.67* 


0.09* 


0 . 00* 


21.00 


5.89* 


3.29* 


1.02* 


0.31* 


0.03* 


0 . 00* 


0.20 


92.84* 


92.32* 


91.27* 


90 . 23* 


88.18* 


82.33* 


1.00 


78.75* 


76.56* 


72.31* 


68.27* 


60.85* 


43. 19* 


2.00 


66.74* 


63.07* 


56.27* 


50.15* 


39.85* 


20.07* 


4.00 


49.55* 


44.25* 


35.22* 


27.98* 


17.66* 


4.48* 


8.00 


28.68* 


22.87* 


14.49* 


9.14* 


3.65* 


0.23* 


12.00 


17. 17* 


12.22* 


6.16* 


3.09* 


0.78* 


0 . 01* 


18.00 


8.17* 


4.91* 


1.76* 


0.62* 


0.08* 


0.00* 


21.00 


5.71* 


3. 15* 


0.95* 


0.28* 


0.03* 


0 .00* 



0 



: 5 



= 10 



= 15 



= 20 
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TABLE A . 2 . b 



Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 



TRANSMITTANCE IN 3.5 UM BAND - MIDLATITUDE SUMMER PROFILE 

*** EMPIRICAL *** 



Range ! Visibility in Kilometers 

in Km ! 40 30 20 15 10 



0 


0.2 


0 


1.0 


0 


2.0 


0 


4.0 


0 


8.0 


0 


12.0 


0 


18.0 


0 


21.0 


5 


0.2 


5 


1.0 


5 


2.0 


5 


4.0 


5 


8.0 


5 


12.0 


5 


18.0 


5 


21.0 


10 


0.2 


10 


1.0 


10 


2.0 


10 


4.0 


10 


8.0 


10 


12.0 


10 


18.0 


10 


21.0 


15 


0.2 


15 


1.0 


15 


2.0 


15 


4.0 


15 


8.0 


15 


12.0 


15 


18.0 


15 


21.0 


20 


0.2 


20 


1.0 


20 


2.0 


20 


4.0 


20 


8. 0 


20 


12.0 


20 


18.0 


20 


21.0 



86.48% 

82.52% 

77.82% 

69.21% 

54.74% 

43.30% 

30.46% 

25.55% 

85.70% 

78.85% 

71.06% 

57.70% 

38.05% 

25.10% 

13.44% 

9.84% 

85.49% 

77.90% 

69.36% 

54.97% 

34.54% 

21.70% 

10.81% 

7.63% 

85.44% 

77.66% 

68.92% 

54.29% 

33.68% 

20.90% 

10 . 21 % 

7.14% 

85.42% 

77.59% 

68.81% 

54.11% 

33.46% 

20.70% 

10.06% 

7.02% 



86. 39% 
82.07% 
76.99% 
67.73% 
52.43% 
40.59% 
27.64% 
22.81% 

85.07% 
76.01% 
66.03% 
49.83% 
28. 38% 
16. 16% 
6.95% 
4.55% 

84.79% 

74.78% 

63.90% 

46.67% 

24.90% 

13.28% 

5.17% 

3.23% 

84.72% 

74.46% 

63.36% 

45.89% 

24.06% 

12.62% 

4.79% 

2.95% 

84.70% 

74.38% 

63.23% 

45.68% 

23.85% 

12.45% 

4.70% 

2.89% 



86.35% 

81.88% 

76.62% 

67.09% 

51.45% 

39.45% 

26.49% 

21.70% 

84.11% 
71.83% 
58.97% 
39.74% 
18. 05% 
8 . 20 % 
2.51% 
1.39% 

83.73% 
70.20% 
56.32% 
36. 25% 
15.02% 
6 . 22 % 
1 . 66 % 
0 . 86 % 

83.63% 

69.78% 

55.65% 

35.40% 

14.32% 

5.79% 

1.49% 

0.76% 

83.60% 
69.68% 
55.48% 
35. 18% 
14.15% 
5.69% 
1.45% 
0.73% 



86. 26% 
81.47% 
75.86% 
65.77% 
49.44% 
37.16% 
24.22% 
19.55% 

83.02% 

67.29% 

51.74% 

30.60% 

10.70% 

3.74% 

0.77% 

0.35% 

82.52% 

65.26% 

48.67% 

27.08% 

8.38% 

2.59% 

0.45% 

0.19% 

82.39% 
64.75% 
47.91% 
26.23% 
7.86% 
2.36% 
0.39% 
0. 16% 

82.35% 

64.61% 

47.71% 

26.02% 

7.74% 

2.30% 

0.37% 

0.15% 



86.09% 

80.65% 

74.34% 

63.16% 

45.60% 

32.92% 

20.19% 

15.81% 

80.69% 

58.34% 

38.89% 

17.29% 

3.42% 

0.67% 

0.06% 

0 . 02 % 

79.93% 
55 . 65% 
35.39% 
14.32% 
2.34% 
0.38% 
0.03% 
0 . 01 % 

79.73% 

54.97% 

34.54% 

13.63% 

2 . 12 % 

0.33% 

0 . 02 % 

0 . 01 % 

79.68% 
54.80% 
34.32% 
13.46% 
2.07% 
0.32% 
0 . 02 % 
0 . 00 % 



5 



85.79% 

79.28% 

71.83% 

58.96% 

39.73% 

26.77% 

14.81% 

11 . 01 % 

74.23% 
38.46% 
16. 90% 
3.26% 
0 . 12 % 
0 . 00 % 
0 . 00 % 
0 . 00 % 

72.83% 
34.95% 
13.96% 
2.23% 
0.06% 
0 . 00 % 
0 . 00 % 
0 . 00 % 

72.47% 

34.09% 

13.28% 

2 . 02 % 

0.05% 

0 . 00 % 

0 . 00 % 

0 . 00 % 

72.37% 

33.87% 

13.11% 

1.96% 

0.04% 

0 . 00 % 

0 . 00 % 

0 . 00 % 
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TAB LB A. 3. a 



TRANSMITTANCE IN 3.5 UM BAND - MIDLATITUDE WINTER PROFILE 

*** LOWTRAN DATA *** 



Range ! 



Visibility in Kilometers 



Wind=0 



Wind~5 



Wind= 10 



Wind=15 



Wind=20 



in Ki 1 


40 


30 


20 


15 


10 


5 


0.2 


97.55* 


97.52* 


97.45* 


97.39* 


97.26* 


96.86* 


1.00 


92.46* 


92.30* 


91.99* 


91 . 68* 


91.06* 


89.22* 


2.00 


88.34* 


88.04* 


87.45* 


86.85* 


85.68* 


82.26* 


4.00 


82.24* 


81.69* 


80.59* 


79.50* 


77.36* 


71.32* 


8.00 


73.25* 


72.28* 


70.35* 


68.45* 


64.82* 


55.09* 


12.00 


66.20* 


64.88* 


62.30* 


59.80* 


55.10* 


43.18* 


18.00 


57.51* 


55.80* 


52.51* 


49.38* 


43.68* 


30.29* 


21.00 


53.76* 


51.90* 


48.34* 


45.00* 


39.00* 


25.45* 


0.2 


96.50* 


96.07* 


95.21* 


94.35* 


92. 66* 


87.78* 


1.00 


87.57* 


85.64* 


81.88* 


78.25* 


71.47* 


54.55* 


2.00 


79.25* 


75.79* 


69.27* 


63.27* 


52.78* 


30.75* 


4.00 


66. 19* 


60.54* 


50.58* 


42. 18* 


29.36* 


9.96* 


8.00 


47.45* 


39.69* 


27.70* 


19.27* 


9.34* 


1.08* 


12.00 


34.51* 


26.40* 


15.40* 


8.93* 


3.01* 


0.12* 


18.00 


21 . 65* 


14.49* 


6.45* 


2.85* 


0.56* 


0.00* 


21.00 


17.19* 


10.76* 


4.19* 


1.62* 


0.24* 


0.00* 


0.2 


96. 28* 


95.77* 


94.75* 


93.73* 


91.72* 


85.99* 


1.00 


86.58* 


84.31* 


79.90* 


75.69* 


67.94* 


49.21* 


2.00 


77.47* 


73.45* 


65 . 98* 


59.21* 


47.69* 


25 . 02* 


4.00 


63.24* 


56.86* 


45.87* 


36.94* 


23.97* 


6.60* 


8.00 


43.32* 


35.01* 


22.79* 


14.78* 


6.22* 


0.47* 


12.00 


30. 11* 


21 . 87* 


11.49* 


6.00* 


1.64* 


0.03* 


18.00 


17.64* 


10.92* 


4. 16* 


1.57* 


0.22* 


0.00* 


21.00 


13.54* 


7.74* 


2.51* 


0.81* 


0.08* 


0 . 00* 


0.2 


96.22* 


95 . 69* 


94. 62* 


93.56* 


91.47* 


85.51* 


1.00 


86.31* 


83.95* 


79.37* 


75.01* 


67.00* 


47 . 83* 


2.00 


76.98* 


72.82* 


65.10* 


58.13* 


46.38* 


23.64* 


4.00 


62.46* 


55.88* 


44.66* 


35. 62* 


22.67* 


5.89* 


8.00 


42.25* 


33. 82* 


21.60* 


13.74* 


5.57* 


0.38* 


12.00 


29.00* 


20 . 77* 


10.60* 


5.38* 


1 . 39* 


0.02* 


18.00 


16.68* 


10.11* 


3.69* 


1.33* 


0. 17* 


0.00* 


21.00 


12.68* 


7.07* 


2.18* 


0.66* 


0.06* 


0.00* 


0.2 


96. 19* 


95 . 65* 


94.56* 


93.48* 


91 . 35* 


85 . 28* 


1.00 


86. 19* 


83.78* 


79. 12* 


74.69* 


66.56* 


47.20* 


2.00 


76.76* 


72.52* 


64.69* 


57.64* 


45.78* 


23 . 02* 


4.00 


62.09* 


55.43* 


44.10* 


35.01* 


22.09* 


5.59* 


8.00 


41.76* 


33.28* 


21.07* 


13.28* 


5.28* 


0.34* 


12.00 


28.49* 


20.27* 


10.21* 


5.11* 


1 . 28* 


0.02* 


18.00 


16.24* 


9.75* 


3.48* 


1.23* 


0.16* 


0.00* 


21.00 


12.30* 


6.78* 


2.04* 


0.61* 


0 . 05* 


0.00* 
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TABLE 4 . 3 . b 



Wind 

Wind 

Wind: 

Wind 

Wind 

Wind 

Wind: 

Wind 

Wind: 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind : 

Wind 

Wind 

Wind : 

Wind 

Wind : 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind: 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 



TRANSMITTANCE IN 3.5 DM BAND - MIDLATITUDE WINTER PROFILE 

*** EMPIRICAL *** 



Range I Visibility in Kiloaeters 

in Kb ! 40 30 20 15 10 



0 


0.2 


0 


1.0 


0 


2.0 


0 


4.0 


0 


8.0 


0 


12.0 


0 


18.0 


0 


21.0 


5 


0.2 


5 


1.0 


5 


2.0 


5 


4.0 


5 


8.0 


5 


12.0 


5 


18.0 


5 


21.0 


10 


0.2 


10 


1.0 


10 


2.0 


10 


4.0 


10 


8.0 


10 


12.0 


10 


18.0 


10 


21.0 


15 


0.2 


15 


1.0 


15 


2.0 


15 


4.0 


15 


8.0 


15 


12.0 


15 


18.0 


15 


21.0 


20 


0.2 


20 


1.0 


20 


2.0 


20 


4.0 


20 


8.0 


20 


12.0 


20 


18.0 


20 


21.0 



93.30% 

91.21% 

88 . 66 % 

83.77% 

74.78% 

66.76% 

56.32% 

51.72% 

92.32% 

86.51% 

79.75% 

67.79% 

48.97% 

35.38% 

21.73% 

17.03% 

92.09% 
85.46% 
77.83% 
64.57% 
44.43% 
30.57% 
17.45% 
13. 18% 

92.04% 

85.22% 

77.40% 

63.85% 

43.45% 

29.57% 

16.60% 

12.43% 

92.03% 

85.17% 

77.30% 

63.69% 

43.23% 

29.34% 

16.41% 

12.27% 



93.24% 
90.93% 
88 . 12 % 
82.75% 
72.98% 
64.37% 
53.31% 
48.52% 

91.65% 

83.44% 

74.19% 

58.67% 

36.68% 

22.94% 

11.34% 

7.97% 

91.33% 

81.98% 

71.63% 

54.68% 

31.87% 

18.57% 

8.26% 

5.51% 

91.26% 

81.65% 

71.06% 

53.81% 

30.86% 

17.70% 

7.69% 

5.06% 

91.24% 

81.58% 

70.92% 

53.61% 

30.63% 

17.50% 

7.56% 

4.97% 



93.23% 

90.87% 

87.99% 

82.52% 

72.57% 

63.83% 

52.64% 

47.81% 

90.61% 

78.80% 

66.18% 

46.67% 

23.22% 

11.55% 

4.05% 

2.40% 

90.14% 
76.77% 
62.81% 
42.04% 
18. 83% 
8.44% 
2.53% 
1.39% 

90.03% 

76.31% 

62.06% 

41.04% 

17.95% 

7.85% 

2.27% 

1 . 22 % 

90.01% 

76.20% 

61.89% 

40.82% 

17.76% 

7.73% 

2 . 22 % 

1.19% 



93.15% 

90.47% 

87.22% 

81.08% 

70.06% 

60.54% 

48.63% 

43.59% 

89.09% 
72.43% 
55.90% 
33.31% 
11.82% 
4.20% 
0.89% 
0 .41% 

88.41% 
69.67% 
51.73% 
28.52% 
8.67% 
2.63% 
0.44% 
0. 18% 

88.25% 

69.05% 

50.82% 

27.53% 

8.07% 

2.37% 

0.38% 

0.15% 

88 . 21 % 
68.91% 
50.61% 
27.30% 
7.94% 
2.31% 
0 . 36% 
0.14% 



92.97% 

89.59% 

85.54% 

77.99% 

64.82% 

53.87% 

40.82% 

35.53% 

85.29% 
58.24% 
36.15% 
13.93% 
2.07% 
0.31% 
0 . 02 % 
0 . 00 % 

84.08% 
54.22% 
31.34% 
10.46% 
1.17% 
0. 13% 
0 . 00 % 
0 . 00 % 

83.81% 
53.34% 
30.33% 
9.80% 
1 . 02 % 
0 . 11 % 
0 . 00 % 
0 . 00 % 

83.75% 
53. 15% 
30.10% 
9.66% 
0.99% 
0 . 10 % 
0 . 00 % 
0 . 00 % 



5 



92.65% 

88.08% 

82.68% 

72.85% 

56.56% 

43.92% 

30.05% 

24.85% 

73.25% 
27.21% 
7.89% 
0 . 66 % 
0 . 00 % 
0 . 00 % 
0 . 00 % 
0 . 00 % 

70.58% 
22.60% 
5.44% 
0.32% 
0 . 00 % 
0 . 00 % 
0 . 00 % 
0 . 00% 

69.99% 
21 . 66 % 
5.00% 
0.27% 
0 . 00% 
0.00% 
0.00% 
0.00% 

69.85% 
21.45% 
4.91% 
0.26% 
0.00% 
0 . 00% 
0. 00% 
0 . 00% 
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TABLE A. 4. a 



TRANSMITTANCE IN 3.5 UM BAND - SUB-ARCTIC SUMMER PROFILE 

*** LOWTRAN DATA *** 



Range ! Visibil 

in Kb I 40 Kb 30 Kb 



Wind=0 0.20 


95 . 52 % 


95.49* 


1.00 


87.64X 


87.49* 


2.00 


81.41* 


81.14* 


4.00 


72.23* 


71.75* 


8.00 


59.00* 


58.22* 


12.00 


49.49* 


48.51* 


18.00 


38.90* 


37.75* 


21.00 


34.70* 


33 . 50* 


Wind=5 0.20 


94.49* 


94.07* 


1.00 


83.00* 


81 . 17* 


2.00 


73.02* 


69.84* 


4.00 


58.11* 


53. 15* 


8.00 


38.19* 


31 . 95* 


12.00 


25.78* 


19.72* 


18.00 


14.62* 


9.79* 


21.00 


11.08* 


6.93* 


Wind=10 0.20 


94.28* 


93.78* 


1.00 


82.08* 


79.93* 


2.00 


71.42* 


67.72* 


4.00 


55 . 58* 


49.98* 


8.00 


34.94* 


28.26* 


12.00 


22.55* 


16.40* 


18.00 


11.97* 


7.42* 


21.00 


8.77* 


5.02* 


Wind= 15 0.20 


94.22* 


93.70* 


1 . 00 


81.83* 


79.60* 


2.00 


70.98* 


67.16* 


4.00 


54.91* 


49.15* 


8.00 


34.10* 


27.32* 


12.00 


21.74* 


15.59* 


18.00 


11.33* 


6.88* 


21.00 


8.23* 


4.60* 


Wind=20 0.20 


94.20* 


93.67* 


1.00 


81.72* 


79.44* 


2.00 


70.78* 


66.89* 


4.00 


54.60* 


48.76* 


8.00 


33.71* 


26.89* 


12.00 


21 . 37* 


15.23* 


18.00 


11.04* 


6.64* 


21.00 


7.98* 


4.41* 



ity in Kiloaeters 



0 

o 

eg 


15 Kb 


10 Kb 


5 Kb 


95.43* 


95.36* 


95.23* 


94.85* 


87.20* 


86.91* 


86.32* 


84.60* 


80.60* 


80.05* 


78.98* 


75.86* 


70.79* 


69.84* 


67.98* 


62.72* 


56. 68* 


55.17* 


52.26* 


44.48* 


46.60* 


44.74* 


41.26* 


32.40* 


35.54* 


33.44* 


29.61* 


20.60* 


31.23* 


29 . 09* 


25.24* 


16.53* 


93.23* 


92.38* 


90.73* 


85 . 95* 


77.60* 


74.16* 


67.73* 


51.69* 


63.83* 


58.29* 


48.63* 


28.32* 


44.40* 


37.03* 


25.77* 


8.74* 


22.30* 


15.51* 


7.51* 


0 . 86* 


11.50* 


6.67* 


2.25* 


0 . 09* 


4.36* 


1.92* 


0.38* 


0.00* 


2.70* 


1.04* 


0.16* 


0 . 00* 


92.78* 


91.79* 


89.83* 


84.23* 


75.77* 


71.79* 


64.45* 


46.73* 


60.85* 


54.63* 


44.03* 


23.15* 


40.36* 


32.52* 


21.13* 


5.84* 


18.42* 


11.96* 


5.05* 


0.39* 


8.63* 


4.52* 


1.24* 


0.03* 


2.83* 


1.07* 


0.15* 


0.00* 


1.63* 


0.53* 


0.05* 


0.00* 


92.66* 


91.62* 


89.58* 


83.77* 


75.28* 


71.15* 


63.58* 


45.46* 


60.06* 


53.66* 


42.85* 


21.90* 


39.31* 


31.38* 


20.01* 


5.23* 


17.48* 


11.14* 


4.53* 


0.31* 


7.98* 


4.06* 


1.05* 


0.02* 


2.52* 


0.91* 


0.12* 


0.00* 


1.42* 


0.44* 


0.04* 


0.00* 


92.60* 


91.55* 


89.47* 


83.55* 


75.05* 


70.86* 


63.18* 


44.87* 


59.70* 


53.22* 


42.31* 


21.34* 


38.84* 


30.86* 


19.51* 


4.96* 


17.06* 


10.77* 


4.30* 


0.28* 


7.69* 


3.86* 


0.97* 


0 .02* 


2.38* 


0.85* 


0.11* 


0.00* 


1.34* 


0.40* 


0.04* 


0.00* 
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TABLE A. 4 . b 



Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind: 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 

Wind 



TRANSMITTANCE IN 3.5 UM BAND - SUB-ARCTIC SUMMER PROFILE 

*** EMPIRICAL *** 



Range ! Visibi 

in Kb ! 40 30 



0 


0.2 


92 . 92% 


92.87% 


0 


1.0 


89.38% 


89.11% 


0 


2.0 


85. 13% 


84.63% 


0 


4.0 


77.24% 


76.33% 


0 


8.0 


63 . 59% 


62.09% 


0 


12.0 


52.35% 


50.51% 


0 


18.0 


39.10% 


37.06% 


0 


21.0 


33.79% 


31.75% 


5 


0.2 


91.95% 


91.33% 


5 


1.0 


84.81% 


81.96% 


5 


2.0 


76.65% 


71.59% 


5 


4.0 


62.62% 


54.63% 


5 


8.0 


41.79% 


31.80% 


5 


12.0 


27.89% 


18.52% 


5 


18.0 


15.20% 


8.23% 


5 


21.0 


11.23% 


5.48% 


10 


0.2 


91.73% 


91.04% 


10 


1.0 


83 . 80% 


80.67% 


10 


2.0 


74.85% 


69.35% 


10 


4.0 


59.70% 


51.25% 


10 


8. 0 


37.99% 


28.00% 


10 


12.0 


24.17% 


15.29% 


10 


18.0 


12.27% 


6.17% 


10 


21.0 


8.74% 


3.92% 


15 


0.2 


91.68% 


90.97% 


15 


1.0 


83.58% 


80.38% 


15 


2.0 


74.45% 


68.86% 


15 


4.0 


59.07% 


50.53% 


15 


8.0 


37.19% 


27.21% 


15 


12.0 


23.41% 


14.66% 


15 


18.0 


11.70% 


5.79% 


15 


21.0 


8.27% 


3.64% 


20 


0.2 


91.67% 


90.96% 


20 


1.0 


83.53% 


80. 32% 


20 


2.0 


74.36% 


68.75% 


20 


4.0 


58.93% 


50.37% 


20 


8.0 


37.01% 


27.04% 


20 


12.0 


23.25% 


14.52% 


20 


18.0 


11.57% 


5.71% 


20 


21.0 


8. 16% 


3.58% 



ity in Kiloaeters 



20 


15 


10 


5 


92.85% 


92.77% 


92.59% 


92.28% 


89.05% 


88.66% 


87.81% 


86.34% 


84.51% 


83.77% 


82.17% 


79.44% 


76.11% 


74.79% 


71.96% 


67.26% 


61.74% 


59.62% 


55.19% 


48.22% 


50.08% 


47.52% 


42.32% 


34.57% 


36.59% 


33.82% 


28.43% 


20.98% 


31.28% 


28.53% 


23.30% 


16.35% 


90.33% 


89.08% 


86.25% 


78. 02% 


77.59% 


72.38% 


61.59% 


37.30% 


64.16% 


55.83% 


40.43% 


14.83% 


43.88% 


33.22% 


17.42% 


2.34% 


20.52% 


11.76% 


3.23% 


0.06% 


9.60% 


4.17% 


0 . 60% 


0.00% 


3.07% 


0.88% 


0.05% 


0.00% 


1.74% 


0.40% 


0.01% 


0.00% 


89. 93% 


88.54% 


85.40% 


76.34% 


75.87% 


70.20% 


58.61% 


33.46% 


61.36% 


52.52% 


36.61% 


11.93% 


40. 12% 


29.40% 


14.29% 


1.52% 


17.15% 


9.21% 


2.17% 


0.02% 


7.34% 


2.89% 


0.33% 


0.00% 


2.05% 


0.51% 


0.02% 


0.00% 


1.08% 


0.21% 


0.00% 


0.00% 


89.84% 


88.42% 


85.21% 


75.97% 


75.50% 


69.72% 


57.97% 


32.66% 


60.75% 


51.81% 


35.81% 


11.37% 


39.33% 


28.61% 


13.67% 


1.38% 


16.48% 


8.72% 


1.99% 


0.02% 


6.91% 


2.66% 


0 . 29% 


0 . 00% 


1.87% 


0.45% 


0.02% 


0.00% 


0.98% 


0.18% 


0.00% 


0.00% 


89. 82% 


88.39% 


85.17% 


75.89% 


75.41% 


69.62% 


57.82% 


32.48% 


60.61% 


51.65% 


35.63% 


11.24% 


39.15% 


28.43% 


13.53% 


1.35% 


16.34% 


8.62% 


1.95% 


0.02% 


6.82% 


2.61% 


0.28% 


0 . 00% 


1.84% 


0.44% 


0.02% 


0.00% 


0.95% 


0.18% 


0.00% 


0.00% 



91 



TABLE A. 5. a 



TRANSMITTANCE IN 3.5 UM BAND - SUB-ARCTIC WINTER 

*** LOWTRAN DATA *** 



Wind=0 



Wind=5 



Wind= 10 



Wind=15 



Wind=20 



Range ! 
in Kb 1 


40 


Visib 

30 


0.20 


98.73* 


98.70* 


1.00 


95.78* 


95.62* 


2.00 


93.00* 


92.68* 


4.00 


88.76* 


88. 16* 


8.00 


82.56* 


81.44* 


12.00 


77.53* 


75.96* 


18.00 


71.13* 


68. 98* 


21.00 


68.30* 


65.90* 


0.20 


97.67* 


97.24* 


1.00 


90.73* 


88.73* 


2.00 


83.45* 


79.81* 


4.00 


71.47* 


65.38* 


8.00 


53.53* 


44.78* 


12.00 


40.47* 


30.98* 


18.00 


26.83* 


17.96* 


21.00 


21.90* 


13.71* 


0.20 


97.44* 


96.92* 


1.00 


89.67* 


87.30* 


2.00 


81.50* 


77.26* 


4.00 


68. 18* 


61.25* 


8. 00 


48.70* 


39. 32* 


12.00 


35.13* 


25.48* 


18.00 


21.70* 


13.40* 


21.00 


17.09* 


9.74* 


0.20 


97.38* 


96.83* 


1 . 00 


89.37* 


86.90* 


2.00 


80.97* 


76.56* 


4.00 


67.29* 


60. 16* 


8.00 


47.44* 


37.92* 


12.00 


33.78* 


24.13* 


18. 00 


20.45* 


12.35* 


21.00 


15.96* 


8.86* 


0.20 


97.35* 


96.79* 


1.00 


89.24* 


86.72* 


2.00 


80.72* 


76.24* 


4.00 


66.88* 


59.65* 


8.00 


46.86* 


37.28* 


12.00 


33. 16* 


23.53* 


18.00 


19.90* 


11.89* 


21.00 


15.45* 


8.47* 



ility in Kiloaeters 



20 


15 


10 


98.63* 


98.56* 


98.43* 


95.29* 


94.96* 


94.30* 


92.05* 


91.41* 


90.15* 


86.96* 


85.76* 


83.42* 


79.23* 


77.06* 


72.91* 


72.89* 


69.92* 


64.34* 


64.84* 


60.92* 


53.78* 


61.31* 


57.01* 


49.29* 


96. 37* 


95.50* 


93.78* 


84.83* 


81.08* 


74.06* 


72.96* 


66.64* 


55 . 60* 


54.63* 


45.57* 


31.73* 


31.27* 


21.76* 


10.55* 


18.07* 


10.49* 


3.54* 


8.00* 


3.54* 


0.69* 


5.34* 


2.06* 


0.31* 


95.88* 


94.84* 


92.80* 


82.71* 


78.33* 


70.26* 


69.36* 


62.20* 


50.05* 


49.37* 


39.71* 


25.71* 


25.54* 


16.52* 


6.92* 


13.34* 


6.94* 


1.88* 


5.08* 


1.91* 


0.27* 


3.14* 


1.00* 


0.10* 


95.74* 


94.66* 


92.53* 


82. 13* 


77.59* 


69.24* 


68.39* 


61.02* 


48.61* 


48.00* 


38.22* 


24.25* 


24.14* 


15 . 30* 


6.16* 


12.26* 


6. 19* 


1.58* 


4.47* 


1.60* 


0.21* 


2.71* 


0.82* 


0.08* 


95.68* 


94.58* 


92.40* 


81 . 86* 


77.24* 


68.77* 


67.94* 


60.48* 


47.94* 


47.37* 


37.54* 


23.59* 


23.51* 


14.76* 


5.83* 


11.78* 


5.86* 


1.46* 


4.21* 


1.48* 


0. 18* 


2.53* 


0.74* 


0.07* 



PROFILE 



5 



98 . 02% 
92 . 38% 
86.51X 
76. 80X 
61.81* 
50 . 22* 
37.08* 
31.94* 

88.85* 
56.54* 
32.41* 
10.78* 
1 . 22 * 
0.14* 
0 . 01 * 
0 . 00 * 

86.97* 
50.81* 
26. 17* 
7.03* 
0.52* 
0.04* 
0 . 00 * 
0 . 00 * 

86.46* 

49.32* 

24.66* 

6.24* 

0.41* 

0.03* 

0 . 00 * 

0 . 00 * 

86 . 22 * 
48.64* 
23.98* 
5.90* 
0.36* 
0 . 02 * 
0 . 00 * 
0 . 00 * 
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TABLE A . 5 . b 



TRANSMITTANCE IN 3.5 UM BAND - SUBARCTIC WINTER PROFILE 

*** EMPIRICAL *** 

Range ! Visibility in Kilometers 



in Km 


40 


30 


20 


15 


10 


5 


Wind= 0 


0.2 


96 . 02% 


95.96* 


95.95* 


95 . 86* 


95.68* 


95.34* 


Wind= 0 


1.0 


94.69* 


94.40* 


94.33* 


93.91* 


92.98* 


91.39* 


Wind= 0 


2.0 


93.05* 


92.47* 


92.34* 


91.52* 


89.73* 


86.67* 


Wind= 0 


4.0 


89.85* 


88.74* 


88.49* 


86.92* 


83.55* 


77.96* 


Wind= 0 


8.0 


83.78* 


81.73* 


81.26* 


78.41* 


72.44* 


63.08* 


Wind= 0 


12.0 


78. 12* 


75.27* 


74.63* 


70.73* 


62.81* 


51.03* 


Wind= 0 


18. 0 


70.34* 


66.52* 


65.67* 


60.59* 


50.71* 


37.14* 


Wind= 0 


21.0 


66.74* 


62.54* 


61.61* 


56.09* 


45.57* 


31.68* 


Wind= 5 


0.2 


95.02* 


94.28* 


93.17* 


91.28* 


86. 18* 


69.38* 


Wind= 5 


1.0 


89.86* 


86.39* 


81.42* 


73.51* 


55.13* 


18.65* 


Wind= 5 


2.0 


83.80* 


77.46* 


68.80* 


56.08* 


31.54* 


3.61* 


Wind = 5 


4.0 


72.87* 


62.26* 


49.12* 


32.64* 


10.32* 


0.14* 


Wind= 5 


8.0 


55. 11* 


40.23* 


25.04* 


11.06* 


1.11* 


0.00* 


Wind= 5 


12.0 


41.68* 


26.00* 


12.76* 


3.75* 


0.12* 


0.00* 


Wind= 5 


18.0 


27.41* 


13.50* 


4.65* 


0.74* 


0 . 00* 


0.00* 


Wind= 5 


21.0 


22.23* 


9.73* 


2.80* 


0.33* 


0 . 00* 


0.00* 


Wind= 10 


0.2 


94.78* 


93.90* 


92.60* 


90. 38* 


84.44* 


65.34* 


Wind=10 


1.0 


88.72* 


84.69* 


78.95* 


69.96* 


49.78* 


13.81* 


Wind=10 


2.0 


81.68* 


74.43* 


64.69* 


50.80* 


25.71* 


1.98* 


Wind=10 


4.0 


69.24* 


57.49* 


43.43* 


26.78* 


6.86* 


0.04* 


Wind=10 


8.0 


49.76* 


34.30* 


19.57* 


7.44* 


0.49* 


0.00* 


Wind= 10 


12.0 


35.76* 


20.46* 


8.82* 


2.07* 


0 . 03* 


0 . 00* 


Wind=10 


18.0 


21.78* 


9.43* 


2.67* 


0.30* 


0 . 00* 


0.00* 


Wind= 10 


21.0 


17.00* 


6.40* 


1.47* 


0. 12* 


0.00* 


0.00* 


Wind= 15 


0.2 


94.72* 


93.81* 


92.46* 


90.17* 


84.02* 


64.39* 


Wind=15 


1.0 


88.44* 


84.28* 


78.36* 


69. 13* 


48.56* 


12.84* 


Wind=15 


2.0 


81.18* 


73.71* 


63.73* 


49.59* 


24.47* 


1.71* 


Wind=15 


4.0 


68. 39* 


56.38* 


42.14* 


25.52* 


6.21* 


0.03* 


Wind=15 


8.0 


48.54* 


32.99* 


18.43* 


6.76* 


0.40* 


0.00* 


Wind=15 


12.0 


34.45* 


19.30* 


8.06* 


1.79* 


0.03* 


0.00* 


Wind=15 


18.0 


20.60* 


8.64* 


2.33* 


0.24* 


0.00* 


0.00* 


Wind=15 


21.0 


15.93* 


5.78* 


1.25* 


0 . 09* 


0 . 00* 


0.00* 


Wind=20 


0.2 


94.71* 


93.79* 


92.42* 


90. 11* 


83.91* 


64. 16* 


Wind=20 


1.0 


88.38* 


84. 18* 


78.22* 


68.92* 


48.26* 


12.61* 


Wind=20 


2.0 


81.06* 


73.53* 


63.49* 


49.30* 


24.17* 


1.65* 


Wind=20 


4.0 


68. 18* 


56.11* 


41.84* 


25.22* 


6.06* 


0.03* 


Wind=20 


8.0 


48.24* 


32.68* 


18. 17* 


6.60* 


0.38* 


0.00* 


Wind=20 


12.0 


34.14* 


19.03* 


7.89* 


1.73* 


0.02* 


0.00* 


Wind=20 


18.0 


20.32* 


8.46* 


2.26* 


0 . 23* 


0 . 00* 


0.00* 


Wind=20 


21.0 


15.67* 


5.64* 


1.21* 


0.08* 


0 . 00* 


0.00* 



93 



TABLE A. 6. a 



TRANSMITTANCE IN 3.5 UM BAND - STANDARD ATMOSPHERE PROFILE 

*** LOWTRAN DATA *** 

Range ! Visibility in Kilometers 

in Km ! 40 30 20 15 10 5 



Wind=0 0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

Wind=5 0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

Wind=10 0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

Wind= 15 0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 

Wind=20 0.20 

1.00 

2.00 

4.00 

8.00 

12.00 

18.00 

21.00 



96 . 67% 


96.64% 


90. 27% 


90. 13% 


85. 25% 


84.99% 


77.78 % 


77.32% 


67.02% 


66.22% 


58.64% 


57.60% 


49.00% 


47.70% 


45.01% 


43.62% 


95.61% 


95.18% 


85.43% 


83.53% 


76.35% 


73 . 00% 


62.39% 


57.04% 


43.12% 


36,03% 


30.26% 


23. 12% 


18.17% 


12.13% 


14. 14% 


8.83% 


95.45% 


94.96% 


84.72% 


82.57% 


75.08% 


71.33% 


60.34% 


54.46% 


40.33% 


32 . 86% 


27.37% 


20.13% 


15.63% 


9.86% 


11.87% 


6.93% 


95.40% 


94.90% 


84.53% 


82.33% 


74.76% 


70.91% 


59.82% 


53 . 82% 


39.64% 


32 . 08% 


26.67% 


19.42% 


15.03% 


9.34% 


11.34% 


6.51% 


95.39% 


94.88% 


84.45% 


82.21% 


74.61% 


70.71% 


59.58% 


53.52% 


39.32% 


31.73% 


26.36% 


19.10% 


14.77% 


9.11% 


11.10% 


6. 32% 



96.58% 


96.52% 


89.86% 


89.59% 


84.48% 


83.97% 


76.40% 


75.47% 


64.65% 


63.09% 


55.55% 


53.57% 


45.19% 


42.78% 


40.95% 


38.42% 


94.32% 


93.46% 


79.84% 


76.28% 


66.68% 


60.87% 


47.60% 


39.66% 


25.10% 


17.42% 


13.44% 


7.77% 


5.37% 


2.36% 


3.42% 


1.31% 


93.98% 


93.01% 


78.42% 


74.44% 


64.34% 


57.97% 


44.30% 


35.97% 


21.74% 


14.33% 


10.84% 


5.80% 


3.89% 


1.52% 


2.34% 


0.79% 


93.90% 


92.89% 


78.06% 


73.97% 


63.74% 


57.24% 


43.49% 


35.07% 


20.95% 


13.63% 


10.25% 


5.38% 


3.58% 


1.36% 


2. 13% 


0.69% 


93.86% 


92.84% 


77.89% 


73.76% 


63.47% 


56.91% 


43. 12% 


34.67% 


20.59% 


13.31% 


9.99% 


5.19% 


3.45% 


1.29% 


2.03% 


0.65% 



96.40% 


96.05% 


89.05% 


87.45% 


82.96% 


80.01% 


73.67% 


68.52% 


60.11% 


52.01% 


49.81% 


40.09% 


38.36% 


27.70% 


33.82% 


23.13% 


91.78% 


86.92% 


69.63% 


53.06% 


50.73% 


24.95% 


27.54% 


9.29% 


8.40% 


0 . 96% 


2.60% 


0.10% 


0.46% 


0 . 00% 


0.19% 


0 . 00% 


91.09% 


85.61% 


67.09% 


49.19% 


47.08% 


25.31% 


23.73% 


6.86% 


6.24% 


0.52% 


1.66% 


0.04% 


0.23% 


0.00% 


0.09% 


0.00% 


90 . 92% 


85 . 28% 


66.44% 


48.23% 


46.18% 


24.34% 


22.83% 


6.34% 


5.77% 


0.45% 


1.48% 


0.03% 


0.20% 


0.00% 


0.07% 


0.00% 


90.84% 


85.12% 


66. 15% 


47 . 80% 


45.77% 


23.90% 


22.43% 


6.12% 


5.57% 


0.41% 


1.41% 


0 . 03% 


0.18% 


0.00% 


0.07% 


0 . 00% 
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TABLE A. 6 . b 



TRANSMITTANCE IN 3.5 UM BAND - STANDARD ATMOSPHERE PROFILE 

*** EMPIRICAL *** 

Range ! Visibility in Kilometers 

in Kh ! 40 30 20 15 10 5 



Wind= 


0 


0.2 


91.41% 


91.36% 


Wind= 


0 


1.0 


88.78% 


88.53% 


Wind= 


0 


2.0 


85.59% 


85.13% 


Wind = 


0 


4.0 


79.56% 


78.70% 


Wind= 


0 


8.0 


68.73% 


67.26% 


Wind= 


0 


12.0 


59.38% 


57.48% 


Wind= 


0 


18.0 


47.69% 


45.42% 


Wind= 


0 


21.0 


42.74% 


40.37% 


Wind= 


5 


0.2 


90.45% 


89.84% 


Wind= 


5 


1.0 


84.21% 


81.41% 


Wind= 


5 


2.0 


77.01% 


71.98% 


Wind= 


5 


4.0 


64.41% 


56.26% 


Wind= 


5 


8.0 


45.06% 


34.38% 


Wind= 


5 


12.0 


31.52% 


21.01% 


Wind= 


5 


18.0 


18.44% 


10.03% 


Wind= 


5 


21.0 


14.11% 


6.93% 


Wind= 


10 


0.2 


90.27% 


89.60% 


Wind = 


10 


1.0 


83.39% 


80.33% 


Wind= 


10 


2.0 


75.53% 


70.07% 


Wind= 


10 


4.0 


61.95% 


53.32% 


Wind= 


10 


8.0 


41.68% 


30.88% 


Wind= 


10 


12.0 


28.04% 


17.88% 


Wind= 


10 


18.0 


15.48% 


7.88% 


Wind= 


10 


21.0 


11.50% 


5.23% 


Wind= 


15 


0.2 


90.24% 


89.55% 


Wind= 


15 


1.0 


83.25% 


80.13% 


Wind= 


15 


2.0 


75.26% 


69.73% 


Wind= 


15 


4.0 


61.51% 


52.80% 


Wind= 


15 


8.0 


41.09% 


30.28% 


Wind = 


15 


12.0 


27.45% 


17.36% 


Wind= 


15 


18.0 


14.99% 


7.54% 


Wind = 


15 


21.0 


11.08% 


4.97% 


Wind= 


20 


0.2 


90.24% 


89.55% 


Wind= 


20 


1.0 


83.22% 


80 . 09% 


Wind= 


20 


2.0 


75.21% 


69.67% 


Wind= 


20 


4.0 


61.43% 


52.71% 


Wind= 


20 


8.0 


40.99% 


30.17% 


Wind= 


20 


12.0 


27.34% 


17.27% 


Wind= 


20 


18.0 


14.90% 


7.48% 


Wind= 


20 


21.0 


11.00% 


4.92% 



91.35% 


91.28% 


91.12% 


90.84% 


88.48% 


88. 13% 


87.37% 


86.06% 


85.02% 


84.36% 


82.91% 


80.43% 


78.50% 


77.28% 


74.65% 


70.26% 


66.92% 


64.86% 


60.52% 


53.61% 


57.06% 


54.44% 


49.06% 


40.90% 


44.91% 


41.86% 


35.81% 


27.26% 


39.85% 


36.70% 


30.60% 


22 . 26% 


88.86% 


87.58% 


84.55% 


75.35% 


77.08% 


71.67% 


60.11% 


33.78% 


64.53% 


55.79% 


39. 23% 


12.39% 


45.22% 


33.80% 


16.72% 


1.67% 


22.20% 


12.40% 


3.04% 


0.03% 


10.90% 


4.55% 


0.55% 


0 . 00% 


3.75% 


1.01% 


0.04% 


0 . 00% 


2.20% 


0.48% 


0.01% 


0 . 00% 


88.52% 


87.10% 


83.77% 


73.72% 


75.60% 


69.74% 


57.37% 


30.28% 


62.07% 


52.82% 


35.75% 


9.96% 


41.84% 


30.30% 


13.88% 


1.08% 


19.01% 


9.97% 


2. 09% 


0.01% 


8.64% 


3.28% 


0.32% 


0 . 00% 


2.65% 


0.62% 


0.02% 


0.00% 


1.46% 


0.27% 


0.00% 


0.00% 


88.46% 


87.01% 


83.63% 


73.42% 


75.33% 


69.39% 


56.89% 


29.68% 


61.63% 


52.29% 


35.15% 


9.57% 


41.25% 


29.70% 


13.42% 


0.99% 


18.48% 


9.58% 


1.95% 


0.01% 


8.28% 


3.09% 


0.28% 


0 . 00% 


2.48% 


0.57% 


0.02% 


0.00% 


1.36% 


0.24% 


0 . 00% 


0 . 00% 


88.45% 


87.00% 


83.60% 


73.37% 


75.28% 


69.33% 


56.80% 


29.57% 


61.55% 


52.20% 


35.04% 


9.50% 


41.15% 


29.59% 


13.33% 


0.98% 


18.39% 


9.51% 


1.93% 


0.01% 


8.22% 


3.06% 


0.28% 


0 . 00% 


2.45% 


0.56% 


0.02% 


0.00% 


1.34% 


0.24% 


0.00% 


0 . 00% 
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TABLE A. 7 



FABLE OF A— COEFF I C I ENTS FOR SIX LOWTRAN PROFILES AT 3.5 urn. 



Tropical Profile 



Visibility 


40 


30 


20 


15 


10 


Wind = 


0 


0 . 85678 1 


0.856471 


0. 856390 


0.856255 


0.856115 


Wi nd = 


5 


0. 856363 


0 . 856405 


0. 856220 


0. 857041 


0.859499 


Wi nd = 


10 


0 • 855773 


0.856587 


0. 856102 


0. 855903 


0. 851515 


Wi nd = 


15 


0 . 856526 


0. 856341 


0. 855761 


0 . 856530 


0.862983 


Wi nd = 


20 


0. 856982 


0.857127 


0. 856171 


0.857050 


0.851565 









Mi dl at i tude 


k Summer Profile 




Vi si bi 1 i ty 


40 


30 


20 


15 


10 


Wi nd = 


0 


0. 875788 


0.875808 


0. 875754 


0. 875834 


0.875812 


Wi nd = 


5 


0. 875747 


0. 875824 


0. 875699 


0. 875477 


0.877011 


Wi nd = 


10 


0.875820 


0. 875716 


0 . 875650 


0. 874656 


0.874052 


Wi nd = 


15 


0. 875753 


0 . 875685 


0 • 875366 


0. 876011 


0. 872420 


Wi nd = 


20 


0. 875777 


0.875885 


0.875886 


0.877431 


0. 859462 









Midlatitude Winter 1 


Prof i 1 e 




Vi si bi 1 i ty 


40 


30 


20 


15 


10 


Wi nd — 


0 


0.938293 


0 . 938300 


0. 938272 


0 . 938302 


0. 938302 


Wi nd = 


5 


0. 938338 


0.938338 


0. 938284 


0 . 937979 


0.938643 


Wind = 


10 


0. 938264 


0.938385 


0.938206 


0.937628 


0.942839 


Wi nd = 


15 


0. 938298 


0. 938291 


0. 938232 


0 . 939094 


0. 942642 


Wi nd = 


20 


0. 938236 


0.938222 


0. 938328 


0.937932 


0. 945244 
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FABLE A. 7 



TABLE OF A— COEFFICIENTS FOR SIX LOWTRAN PROFILES AT 3.5 um. 



Subarctic Summer Profile 



Vi si bi 1 i ty 


40 


30 


20 


15 


10 


Wind = 0 
Wind = 5 
Wind = 10 
Wind = 15 
Wind — 20 


0 . 900942 
0. 900908 
0. 900920 
0.900901 
O. 901048 


0. 900967 
0. 900991 
0.900967 
0. 900910 
0.900953 


0 . 900928 
0 . 900808 
0. 901247 
0.901037 
0.900717 


0.900916 
0.901004 
0.900275 
0. 900193 
0. 900709 


0.900932 
0. 897534 
0. 910749 
0 . 90337 0 
0. 889788 






Subarcti c 


Winter Pro-file 




Vi si bi 1 i ty 


40 


30 


20 


15 


10 


Wind = 0 
Wind = 5 
Wind = 10 
Wind = 15 
Wind = 20 


0.963535 
0.963548 
0. 963594 
0 . 963558 
0. 963587 


0. 963535 
0. 963644 
0.963567 
0. 963605 
0.963634 


0.963587 
0. 963652 
0. 963533 
0 . 963507 
0. 963474 


0.963549 
0.963692 
0. 963565 
0. 963501 
0. 963308 


0. 963600 
0. 963209 
0. 965684 
0. 956474 
0. 956391 






U.S. Standard Atmosphere Pro-file 


Vi si bi 1 1 ty 


40 


30 


20 


15 


10 


Wind = 0 
Wind = 5 
Wind = 10 
Wind = 15 
Wind = 20 


0. 92081 1 
0.920841 
0. 920800 
0. 920802 
0.920836 


0.920821 
0. 920853 
0.920763 
0. 920818 
0. 920866 


0. 920850 
0.920746 
0. 920945 
0.920675 
0. 920934 


0. 920789 
0.920921 
0. 920259 
0.920402 
0. 920379 


0. 920866 
0.922365 
0. 919745 
0. 923306 
0.914907 
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TABLE A. 8 



FABLE OF B— COEFFICIENTS FOR SIX LOWTRAN PROFILES AT 3.5 urn. 



Tropical Profile 



Vi si bi 1 i ty 


40 


30 


20 


1 5 


10 


Wind = 


0 


-0.078300 


-0. 079883 


-0. 083265 


-0.086652 


-0. 093444 


Wi nd = 


5 


-0. 132564 


-0. 154864 


-0. 199849 


-0. 245554 


-0.337061 


Wi nd = 


10 


-0. 143678 


-0. 170376 


-0.223794 


-0. 277764 


-0.384270 


Wi nd = 


15 


-0. 146815 


-0. 174519 


-O. 230264 


-0. 286954 


-0.401796 


Wind = 


20 


-0. 148416 


-0. 176768 


-0. 233486 


-0.291292 


-0. 404436 







Midlatitude Summer 1 


Prof i 1 e 




Visibility 


40 30 


20 


15 


10 


Wi nd = 


0 


-0.064275 -0.065945 


-0. 069313 


-0.07271 1 


-0.079494 


Wi nd = 


5 


-0.118579 -0.140899 


—0. 185840 


-0. 231 112 


-0. 322254 


Wi nd = 


10 


-0.129858 -0.156395 


-0. 209959 


-0.263693 


-0.371488 


Wi nd = 


1 5 


-0. 132919 -0. 160663 


-0. 216479 


-0. 273163 


-0.384626 


Wi nd = 


20 


-0.134394 -0.162714 


-0. 219772 


-0. 277832 


-0.386820 





Midlatitude Winter 1 


Prof lie 




Visibility 


40 30 


20 


15 


10 


Wind = 0 
Wind = 5 
Wind = 10 
Wind = 15 
Wind = 20 


-0.028635 -0.030311 
-0.082925 -0.105237 
-0.094293 -0.120931 
-0 . 0974 16 -0.1 25235 
-0.098874 -0.127235 


-0. 033691 
-0. 150164 
-0. 174555 
-0. 181256 
-0. 184428 


-0.037106 
-0. 195458 
-0.228508 
-0. 238064 
-0. 242007 


-0.043922 
-0.286223 
-0.338264 
-0. 352103 
—0 . 35937 6 
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TABLE A. 8 



TABLE OF B— COEFFICIENTS FOR SIX LOWTRAN PROFILES AT 3.5 um. 



Subarctic Summer Pro-file 



Vi sibi 1 i ty 


40 30 


20 


15 


10 


Wind = 0 
Wind = 5 
Wind = 10 
Wind = 15 
Wind = 20 


-0.048863 -0.050536 
-0. 103225 -0. 125554 
-0.11 4354 -0 . 1 409 1 4 
-0.117388 -0.145095 
-0.118852 -0.147091 


-0. 053881 
-0. 170457 
-0. 194472 
-0. 201005 
-0. 203899 


—0 . 057260 
-0.215911 
-0.248144 
-0. 257037 
-0. 261375 


-0.064021 
-0. 305416 
-0.359279 
-0.370669 
-0.372554 




Subarct i c 


Winter Profile 




Vi si bi 1 i ty 


40 30 


20 


15 


10 


Wind = 0 
Wind = 5 
Wind = 10 
Wind = 15 
Wind = 20 


-0.017768 -0.019474 
-0.071933 -0.094238 
-0.083740 -0.110510 
-0 . 087009 -0.11 5030 
-0.088548 -0.117162 


-0 . 0229 1 2 
-0. 139144 
-0. 164416 
-0. 171452 
-0. 174745 


-0.026375 
-0. 184483 
-0. 218868 
-0. 228422 
-0.233146 


-0. 033305 
-0.274853 
-0.328210 
-0.340038 
-0.346934 




U.S. Standard Atmosphere Profile 


Vi sibi 1 i ty 


40 30 


20 


15 


10 


Wind = 0 
Wind = 5 
Wind = 10 
Wind = 15 
Wind = 20 


-0.036796 -0.038290 
-0.091927 -0.004360 
-0. 100272 -0. 125884 
-0.102441 -0.128875 
-O. 103449 -0. 130276 


-0. 041297 
-0. 159548 
-0. 177567 
-0. 182099 
-0. 184322 


-0.044335 
-0. 205238 
-0.229439 
-0. 235793 
-0. 238687 


-0. 050403 
-0.296848 
-0. 333287 
-0. 344040 
— 0. 345855 
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APPENDIX B 



COMPUTER CODING TO INSTITUTE 
ATMOSPHERIC TRANSMITTANCE CALCULATION 

The -following FORTRAN code is designed to institute the 
calculation of atmospheric transmittance on an incremental 
range movement basis. As the missile closes its target (s) 
periodic recal cul at i on o-f acquisition pr obab i 1 i t i es based on 
updated range and atmospheric tr ansmi t tance are necessary to 
real i stical 1 y model an engagement. Listed here are the 
subroutines necessary to calculate the tr ansmi ttance based 
on the value of the current range value in the model. 
Although these subroutines are written in FORTRAN, they are 
simple and easily convertable to other languages. Comments 
regarding user alternatives are contained therein. 



C First thing to do is allow the user to select which 
C Atmospheric Profile is desired- Then pass this to the 
C ATMOSP subroutine to calculate specific coef f i ci ents . 

C This input may be converted to disk— read input or 

C NAMELIST input if 1 n t er act i veness is not desired. 

INTEGER ATPROF 



111 



READ (111 
FORMAT ( ' 
1 ' 1 . 

2 ' 2 . 

3 3 . 

4 ' 4 . 

5 5 . 

6 ' 6 . 



,*) ATPROF 

Which Atmospheric Profile 
Tropical ' / 

Mi dl ati tude Summer ' / 

Mi dl at i tude Winter'/ 
Subarct i c Summer ' / 

Subarcti c Winter'/ 

U.S. Standard Atmosphere', 



i s 



I 1 ) 



desired 



?'// 



CALL ATMOSP (ATPROF) 



100 



C The values of Wind and Visibility must already have 

C been input to the program prior to the above CALL. 

C ATMOSP calculates the 'A' and *B' coefficients which 

C are necessary to refer to the TRANSh function. 



SUBROUTINE ATMOSP (II) 

DIMENSION A1 (6) ,A2(6) ,A3(6) ,A4(6) ,C1 (6) ,C2(6) , 
1 D1 (6) , D2(6) , SIGMA (6) 

COMMON ACOEFF , BCOEFF , V I S I B , VW I ND 

OPEN ( UN I T=9 , F I LE= ' ATMPROF . DAT ' , STATUS= ' OLD ' ) 

DO 1 1=1,6 



READ (9,*) 


A1 ( I ) 


READ (9 , *) 


A2 ( I ) 


READ (9,*) 


A3( I) 


READ (9, *) 


A4 ( I ) 


READ (9 , *■) 


Cl (I) 


READ (9 , *) 


C2 ( I ) 


READ (9,*) 


D1 ( I ) 


READ (9,*) 


D2 ( I ) 


READ (9 , *) 


SIGMA(I) 



1 CONTINUE 

CLOSE < UNI T=9,STATUS=' KEEP ' ) 

ACOEFF = SIGMA (II) 

VCOEFF = Al(II) + A2(II)*VISIB + A3 ( 1 1 ) *VISIB**2. 

# + A4 (II) *VISIB**3. 

WCOEFF = 1. + D1(II)*(1.— EXP (D2 (II) *VWIND) ) 

IF (VWIND.LE.O. 1) THEN 

VWC0EF=1.0 './* No correlation if 

GO TO 2 './* wind near zero. 

END IF 

VWCOEF = Cl (II) * VISIB ** C2 (II) 

2 BCOEFF = VCOEFF * WCOEFF * VWCOEF 

RETURN 

END 

C The two necessary coefficients (A or 0, B or «) are 
C computed and global for referencing by TRANSM for 
C transmittance calculation based on current range 
C values. Now the function to perform equation 4.3. 
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FUNCTION TRANSM (RNG , AT ) 

C The current range value is passed in as RNG, and AT is 
C the Atmospheric transmittance value passed back out. 

COMMON RO , ACOEFF , BCOEFF 

AT = ACOEFF * EXP < BCOEFF*RNG) 

RETURN 

END 

The Tile 'ATMPROF.DAT' referred to in subroutine ATMOSP 
above is an ASCII data Tile which contains all oT the curve 
Titting coeTTicients listed in table 4.8. It is easily 
constructed with an editor, keeping in mind that iT a name 
change is necessary, it must be perTormed in both places. 
The Tormat is shown here and may be completed by reTerring 
to table 4.8. Ensure that no comments, blank lines or 
leading spaces occur when constructing the Tile. 



-0. 139154 


! /* 


Tropi cal 


A1 


6.2887E-03 


! /* 




A2 


“2. 2354E-04 


! /* 




A3 


2. 609 1 E— 06 


! /* 




A4 


13. 74105 


! /* 




Cl 


-0.716806 


! /* 




C2 


0.895 


! /* 




D1 


-0.30 


1 /* 




D2 


0. 9383 


! /* 




0 


-0. 1249852 


! /* 


Midi at Summer 


A 1 



etc. Tor others in the order speciTied. 
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